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CHAPTER  I 


INTRODUCTION 


This  report  is  concerned  with  the  problem  of  calculating  the  drag 
force  on  an  object  moving  in  a  gaseous  medium.  An  important  parameter 
that  governs  the  properties  of  aerodynamic  flows  in  rarefied  gases  is 
the  Knudsen  number  K.  This  Knudsen  number  is  defined  as  the  ratio  of 
the  mean  free  path  X  of  the  molecules  in  the  gas  and  a  length  L  which 
measures  the  size  of  the  object.  In  order  to  describe  the  aerodynamic 
processes  when  the  Knudsen  number  is  larger  than  unity,  one  needs  to 
consider  methods  based  on  the  kinetic  theory  of  gases.  That  is,  in 
rarefied  gas  dynamics,  the  drag  force  can  be  considered  as  resulting 
from  collisions  of  individual  molecules  with  the  object. 

The  regime  in  which  the  molecular  mean  free  path  is  much  larger  than 
the  size  of  the  object,  i.e.  the  limit  of  infinite  Knudsen  number,  is 
usually  referred  to  as  the  free  molecular  flow  regime.  In  this  regime 
molecules  that  are  reflected  from  the  object  collide  on  the  average  with 
oncoming  molecules  at  large  distances  from  the  object.  Hence  the  process 
is  here  completely  determined  by  the  interaction  of  independent  gas 
molecules  and  the  object,  while  collisions  between  the  molecules  may  be 
neglected.  Aside  from  the  detailed  nature  of  the  interaction  between  gas 
molecules  and  the  object,  the  physics  of  the  free  molecular  flow  regime 
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is  well  understood  and  practical  calculations  abound  in  the  literature. 
Reviews  of  Schaaf  [1]  and  Schaaf  and  Chambre  [2]  give  examples  of  free 
molecular  flow  calculations  together  with  bibliographies. 

In  this  report  we  focus  our  attention  on  the  nearly  free  molecular 
flow  regime  where  the  Knudsen  number  is  greater  than  one,  but  definitely 
finite.  For  such  flows  the  molecular  mean  free  path  is  larger  than  a 
typical  dimension  of  the  object,  but  small  enough  so  that  collisions 
between  the  gas  molecules  cannot  be  neglected. 

The  drag  coefficient  CD  of  an  object  in  a  gas  stream  is  defined  as 


c  =  2. 

D  U„ 


(1-1) 


where  E  is  the  magnitude  of  the  force  exerted  on  the  object  and  U  is  the 

& 

incident  kinetic  energy.  In  the  nearly  free  molecular  flow  regime  K  1«1 
the  drag  coefficient  Cp  may  be  written  in  the  form 


C  =  C  +C.  K  +  .  .  . 
D  o  1 


(1-2) 


where  C  is  the  drag  coefficient  in  the  free  molecular  flow  regime 
o 

K  ■L-»-0  and  the  coefficient  of  a  correction  term  which  is  inversely 
proportional  to  the  Knudsen  number  K. 

t 

The  drag  force  can  be  evaluated  theoretically  by  solving  the 
Boltzmann  equation  subject  to  the  appropriate  boundary  conditions  imposed 
by  the  presence  of  the  object.  The  Boltzmann  equation  is  a  nonlinear 
integro-differential  equation  that  describes  the  rate  of  change  of  the 
single-particle  distribution  function  f  of  a  dilute  gas  [3] 


»f<*  „  3f<?  ;t) 

- it - +V1 - if -  =J(££) 


(1-3) 


“>■  *> 

where  r^,  v^  represent  the  position  and  velocity  of  a  molecule  labeled  1 
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and  where  J  (ff )  is  a  collision  term  containing  a  time  independent  integral 
operator  acting  on  the  function  f (r^ ,v^> t) .  The  term  J(ff)  accounts  for 
the  change  in  f  due  to  binary  collisions  between  the  molecules  and  the 
gas.  An  approximate  solution  of  the  Boltzmann  equation  was  obtained  by 
Liu  et  al.  for  a  sphere  in  a  gas  stream  of  Maxwellian  molecules,  i.e. 
molecules  that  repel  each  other  with  a  force  that  varies  as  the  inverse 
fifth  power  of  the  intermolecular  separation  [4] . 

Since  it  has  proven  to  be  difficult  to  solve  the  full  nonlinear 
Boltzmann  equation,  some  of  the  more  authoritative  studies  are  restricted 
to  flows  at  low  velocities  where  the  linearized  Boltzmann  equation  may  be 
used  [5] .  Many  authors  have  replaced  the  Boltzmann  equation  with  a  model 
equation  proposed  by  Bhatnagar,  Gross  and  Krook  [6]  . 


8_f 

3t 


+  v  •  =  v(f  -f) 

3?1 


(1-4) 


where  f  is  the  local  Maxwell  distribution  and  v  an  adjustable  parameter 
representing  a  velocity  independent  collision  frequency.  Willis  has  solved 
the  BGK  equation  using  a  Knudsen  iteration  method  [7] .  Rose  has  used  a 
Fourier  transform  technique  to  deduce  the  drag  from  the  BGK  equation  [8] . 
Cercignani  and  coworkers  have  evaluated  the  drag  of  a  sphere  from  the  BGK 
equation  using  a  variational  method  [9] .  For  a  review  of  the  use  of  the 
Boltzmann  equation  and  the  BGK  equation  in  rarefied  gas  dynamics  the 
reader  is  referred  to  the  books  of  Kogan  [10]  and  Cercignani  [11] .  Although 
judicious  application  of  the  BGK  equation  has  yielded  some  encouraging 
results,  the  model  equation  does  not  have  any  really  predictive  power  [12]. 
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Another  approach,  encountered  in  the  literature,  is  the  first  col- 
lision  model.  The  first  collision  model  assumes  that  the  drag  on  the 
object  is  the  sum  of  the  free  molecular  flow  drag  Cq  and  a  correction 
term  which  results  from  collisions  between  molecules  reflected  from  the 
surface  of  the  object  and  molecules  in  the  incident  beam.  It  was  used  by 
Lunc  and  Lubonski  [13]  for  calculating  the  drag  of  a  strip,  by  Baker  and 
Charwat  [14]  for  calculating  the  drag  of  a  sphere,  and  by  Perepukhov  for 
calculating  the  drag  of  a  sphere  and  a  cone  [15] .  The  first  collision 
model  is  a  phenomenological  theory  and  it  is  not  derived  from  first 
principles.  Clearly  as  the  Knudsen  number  decreases,  more  and  more  colli¬ 
sions  between  molecules  must  be  considered. 

A  systematic  approach  to  account  for  the  effect  of  molecular  colli¬ 
sions  on  the  drag  was  developed  by  Dorfman  and  coworkers  [16,17].  In 
this  approach  the  object  is  treated  as  a  heavy  particle  and  the  dynamical 
evolution  of  the  system  of  gas  and  object  is  treated  with  the  aid  of 
the  same  techniques  used  earlier  to  derive  the  generalized  Boltzmann 
equation  for  a  moderately  dense  gas  from  the  Liouville  equation.  In 
general  the  theory  leads  to  a  density  expansion  for  the  drag  force  on  the 
object.  In  the  limit  that  the  mean  free  path  of  the  molecules  is  large 
compared  to  the  size  of  the  molecules  themselves,  this  expansion  reduces 
to  an  expansion  in  the  inverse  Knudsen  number.  The  coefficients  in  this 
expansion  are  given  by  integrals  related  to  sequences  of  successive 
collisions  among  the  molecules  and  the  object.  The  same  collision  inte¬ 
grals  can  also  be  obtained  by  solving  the  Boltzmann  evocation  with  a 
modified  Knudsen  number  iteration  procedure  as  shown  by  Kelly  and  Sengers 
[18].  The  derivation  of  the  collision  integrals  for  the  first  inverse 
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Knudsen  number  correction  term  in  (1-2)  is  reviewed  in  Chapter  II.  A 
preliminary  analysis  of  the  collision  sequences  that  enter  into  the 
evaluation  of  the  coefficient  in  (1-2)  was  made  by  Kuperman  [19] . 

In  principle,  our  method  of  calculating  the  drag  coefficient  from 
collision  integrals  can  be  applied  to  objects  of  any  shape  in  a  gas 
stream  of  molecules  with  any  interaction  potential  of  finite  range 
given  the  interaction  mechanism  between  the  molecules  and  the  object. 
However,  in  our  current  studies  of  the  nature  of  the  collision  integrals 
we  have  introduced  the  following  approximations: 

1)  .  The  molecules  that  strike  the  object  do  not  stick  to  it,  but  are 

re-emitted  after  a  time  short  compared  to  the  mean  free  time  of 
the  molecules. 

2)  .  The  molecules  are  re-emitted  diffusively  with  a  temperature  T 

corresponding  to  the  temperature  of  the  object  which  is  assumed 
to  be  the  same  as  the  temperature  of  the  molecules  in  the  gas 
stream.  We  are  thus  neglecting  any  heat  transfer  effects. 

3)  .  The  molecules  in  the  gas  stream  are  assumed  to  interact  as  hard 

spheres  with  mass  m  and  diameter  a . 

The  drag  force  on  the  object  not  only  depends  on  the  Knudsen 
number  of  the  system,  but  also  on  the  Mach  number  of  the  gas  stream. 

The  Mach  number  M  is  defined  as  the  ratio  of  the  flow  velocity  V  relative 
to  the  sound  velocity.  For  the  problem  at  hand,  instead  of  the  Mach 
number,  we  find  it  more  convenient  to  use  as  a  parameter  the  speed  ratio. 
The  speed  ratio  S  is  defined  as  the  ratio  of  the  flow  velocity  V  relative 

1 f 

to  the  thermal  velocity  (2kT^n)  where  k  is  Boltzmann's  constant: 

S  =  V(m/2W)  2.  (1-5) 
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The  speed  ratio  S  is  directly  proportional  to  the  Mach  number  M 

S  =  M(y/2)^  ,  (1-6) 

where  Y  is  the  specific  heat  ratio  Cp/cy  [1]  . 

In  a  preceding  research  effort  Kuperman  analyzed  the  collision 
integrals  in  the  limits  of  zero  and  infinite  Mach  number  [19] .  The 
research  effort  reported  here  is  concerned  with  a  study  of  the  feasibility 
of  evaluating  these  collision  integrals  for  all  Mach  numbers  or  speed 
ratios.  In  particular,  we  shall  present  results  obtained  for  the  drag 
coefficient  as  a  function  of  the  speed  ratio  for  a  disc  and  a  sphere 
in  Chapters  III  and  IV,  respectively.  We  were  informed  at  the  initial 
stage  of  this  research  effort  that  Willis  had  obtained  from  the  Boltzmann 
equation  a  completely  analytic  solution  for  the  drag  coefficient  of  a 
disc  in  the  infinite  Mach  number  limit  [20,21].  Therefore,  by  applying 
our  method  to  a  disc  we  shall  be  able  to  compare  our  method  with  that 
of  Willis  without  any  uncertainty  associated  with  the  finite  precision 
of  any  numerical  quadrature  method  and  to  extend  his  result  to  finite 
Mach  numbers.  The  sphere  drag  was  selected  since  we  had  previously 
studied  the  sphere  drag  in  the  zero  and  infinite  Mach  number  limit  [16,19] 
and  since  other  investigators  have  studied  the  sphere  drag  by  approxi¬ 
mate  methods  [4,7,8,14,15]. 
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CHAPTER  II 


DERIVATION  OF  COLLISION  INTEGRALS  FOR  THE  DRAG  FORCE 
IN  THE  NEARLY  FREE  MOLECULAR  FLOW  REGIME^ 


2.1  Introduction 


A  theory  of  the  force  on  an  object  in  a  gas  stream  in  analogy 
with  current  procedures  in  the  kinetic  theory  of  moderately  dense 
gases  was  proposed  by  Dorfman  and  Sengers  [22]  and  formulated  by 
McClure  and  Dorfman  [17] .  It  is  the  purpose  of  this  chapter  to  review 
the  theory  and  to  derive  the  general  form  of  the  collision  integrals 
for  the  first  inverse  Knudsen  number  correction  to  the  drag  force. 

These  formulas  will  then  be  analyzed  and  evaluated  as  a  function  of  the 
speed  ratio  in  the  subsequent  chapters. 

Let  us  consider  a  system  of  N  gas  molecules  and  a  macroscopic 
object  in  a  volume  fl.  The  object  is  at  rest  and  located  at  the  origin 
of  the  coordinate  system.  We  shall  consider  the  limit  in  which  the  ra¬ 
tio  of  the  mass  of  the  object  and  the  mass  of  the  molecules  becomes  infinite. 
The  momentum  and  position  vectors  of  molecule  i  are  indicated  by  p^  and 
respectively.  The  total  momentum  of  the  gas  at  time  t  is  given  by 
P(t)  =^P±(t)  .  (2-1) 


This  chapter  was  prepared  in  collaboration  with  Professor  J.  R.  Dorfman, 
Dr.  C.  F.  McClure  and  Dr.  W.  A.  Kuperman  [17,19]. 
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In  order  to  evaluate  the  total  momentum  we  use  the  methods  of  non¬ 
equilibrium  statistical  mechanics  with  the  state  of  the  gas  being 

characterized  by  an  N-particle  distribution  function  D  (T  ;X;t) .  Here 

N  N 

_  _  -V  -V 

r  =  (x.  f  •  *.xm,)=  (p.  # r.  #  •  .pMfrw)  represents  the  6N-dimensional  phase  space 
N  1  N  11  N  N 

of  the  gas,  X  denotes  the  macroscopic  object  and  ;X;t)  is  defined 

N  N 

such  that  D„(r  ;x?t)dT  is  the  probability  that  the  phases  of  the  particles 

N  N  N 

in  the  presence  of  the  object  will  lie  between  T  and  T  +  dT  at  time  t. 

N  N  N 

The  distribution  function  is  normalized  such  that 


'  1. 


The  expectation  value  <P(t)>  is  then  determined  by 


(2-2) 


N 

<!(«>  .  KJXW^1  •  (2-3) 

i"=l 

Neglecting  wall  effects  that  disappear  in  the  thermodynamic  limit 
(N-*00,  N/S2=n) ,  the  force  E(t)  on  the  object  is  equal  to  the  negative 

of  the  time  rate  of  change  of  this  total  momentum  [17] 


E(t)  =  -  ££  <  P(t)> 


(2-4) 


The  time  evolution  of  the  N-particle  distribution  function 
DN(rN;t)  in  the  absence  of.  a  foreign  object  is  governed  by  the  Liouville 
equation 


3tDN(IN;t*  +  DN^N;tJ  “  0 


where  LN  is  the  Liouville  operator 


(2-5) 
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vv-v  1 9ii 

i=i  ri  ±<j 

Here  0^  is  a  differential  operator 

S^i-i  3  34 ,  ,  3 

0.  .  =  — ±2  +  -13 J  . - 

13  ??,  »?±  ajj 


(2-6) 


(2-7) 


where  j —  ^(^j)  =  <Mlri“rj|)  is  the  intermolecular  potential  between 
molecules  i  and  j .  The  solution  of  the  Liouville  equation  (2r5)  may  be 
written  formally  as 


dnuv*>  *  •■“"<r,,>vv  °>  • 


(2-8) 


In  the  past  decade  methods  have  been  developed  for  the  evaluation 
of  the  average  value  of  a  phase  function  A(x1...x  ) 

<A>  =  |drN  A(rN)Dh(rN,t)  .  (2-9) 

In  particular,  Cohen  and  Green  have  formulated  a  cluster  expansion  of 
the  streaming  operator  exp(-tLN)  that  leads  to  a  density  expansion  for 
<A>  [23,24].  In  order  to  apply  the  same  method  to  the  drag  problem  we 
need  to  extend  the  concept  of  streaming  operator  or  time-displacement 
operator  to  a  system  of  molecules  in  the  presence  of  a  foreign  object. 


2.2  Time  Evolution  of  a  System  of  Molecules  in  the  Presence  of  an  Object. 

A  description  of  the  time  evolution  of  a  system  of  molecules  in  the 
presence  of  a  foreign  object  must  take  into  account  the  fact  that  the 
interaction  of  the  molecules  and  the  surface  of  the  object  is  usually  of 
a  stochastic  nature.  A  suitable  time-displacement  operator  for  such  a 
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system  was  formulated  by  McClure  and  Dorfman  [17]. 

Let  II  .  .  (T ' ;  t'l  T„;t)  be  the  transition  probability  density  that 
t  — t  N  N 

takes  the  system  from  the  phase  IV  at  time  t  to  the  phase  T '  at  time  t\ 

N  N 


The  average  total  momentum  of  the  system  is  given  by 


<P(t»-  fdT  ?ffN)  DB(rH,x,t) 

i  PQ1^)  =  \  P,^.  In  terms  of  t 


C2-10) 


where  P(T 


rewrite  (2-10)  as 


In  terms  of  the  transition  probability  II  we  may 


<p(t]>  =  JarH  ?<rN>Jar- 

Introducing  a  time  displacement  operator  as 

Ttf<rS>  =/V(IVVIVtlrS;o> 


(2-11)  becomes 


<p(ti>=ydr»  Vr“iX;0)  Tt  >  • 


(2-11) 


(2-12) 


(2-13) 


In  order  to  formulate  the  time-displacement  operator  T  we  need  to 
specify  the  interaction  between  the  molecules  and  the  surface  of  the 
object.  For  this  purpose  we  assume  that  a  gas  molecule  after  colliding 
with  the  object  with  an  incoming  velocity  v,  will  be  re-emitted  from 
the  surface  with  a  velocity  v'  with  a  probability  n(v,|v)dv'.  We 
assume  that  all  molecules  impinging  on  the  object  will  be  re-emitted 


again,  so  that  the  transition  probability  is  normalized  as 


^dv'  n(v’|v)  = 


(2-14) 


To  account  for  the  interaction  of  a  molecule,  labeled  1,  with  the 
surface  of  the  sphere  it  is  convenient  to  introduce  an  operator  T(1X), 
defined  as  [17,19] 


(IX)  =  J* dv'^dA  r|(v'jv)  |  v*nj  6s  (r-R)  (f^-1)  , 


(2-15) 


v.n<o 
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Here  it  is  assumed  that  the  molecule  with  incoming  velocity  v 
strikes  the  surface  of  the  object  at  the  position  R  inside  the  two- 
dimensional  surface  element  dA.  The  symbol  5  (r-R)  represents  a  three- 
dimensional  delta  function.  The  operator  is  a  rotation  operator  that 
transforms  the  velocity  v  before  the  collision  into  the  velocity  v'  after 
the  collision  of  the  molecule  with  the  object.  The  vector  n  is  a  unit 
vector  in  the  direction  of  the  outward  normal  and  the  condition  v*n<o 
indicates  that  the  molecule  strikes  the  surface  from  the  outside.  The 
integrations  extend  over  all  velocities  v1  and  over  the  entire  surface 
of  the  object. 

We  shall  refer  to  T(1X) ,  defined  in  (2-15),  as  the  binary  collision 
operator  for  a  molecule  and  the  object.  It  may  be  considered  as  the  sum  of 
two  operators 

T  (IX)  =  TX(1X)  +  Tn(lX)  ,  (2-16) 

where 


T"*" (IX )  =  /<£  n(v'  |v)  |v*n|53  (r-R)RR  , 


v*n<o 


T  (IX)  =  - J  dA|vn|S3  (r-R) 


v*n<o 


(2 -16a) 


(2 -16b) 


The  meaning  of  the  operators  T^dX)  and  Tn(lX)  is  illustrated  in  Fig.  1. 

In  this  figure  the  circle  indicates  the  surface  of  the  object  and  the  lines 
represent  the  trajectory  of  the  molecule.  The  operators  T1(1X)  and  Tn(lX) 
are  only  different  from  zero,  if  the  molecule  impinges  on  the  object.  The 
operator  T1(1X)  transforms  the  incident  velocity  v^=v  of  the  molecule  into 
the  velocity  v1  of  the  molecule  after  being  re-emitted  from  the  object, 
taking  into  account  the  appropriate  probability  distribution  for  v' .  We 
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refer  to  this  process  as  an  "interacting"  collision  between  the  molecule  and 
the  object  as  shown  in  Fig.  la.  The  operator  Tn(lX)  requires  also  that  the 

41 

molecule  collides  with  the  object,  but  it  does  not  change  the  velocity  v. 

The  operator  Tn(lX),  therefore,  corresponds  to  a  "non- interacting"  collision 
in  which  the  velocity  of  the  molecular  after  the  collision  is  equal  to  the 
velocity  prior  to  the  collision,  as  shown  in  Fig.  lb. 

The  operator  T(1X)  is  a  generalization  of  the  binary  collision  opera¬ 
tors  T(ij)  introduced  by  Ernst  et  al.  to  describe  the  interaction  between 
two  hard  sphere  molecules  [25].  These  binary  collision  operators  T(ij) 


were  defined  as 

T (ij)  a  T3" Ci j )  +  Tn(ij)  , 


(2-17) 


where  T  tij)  corresponds  to  an  interacting  collision  between  two  hard 
sphere  molecules  i  and  j 


T;L(ij)  =  a2  fdcr. .  |v. .  *a. .  |63  (r. 

/  ID1  ID  I]1  ID  ID  cr±j  , 


(2-.17a) 


v. . -a.  ,<o 
ID  ij 


J*,, 


and  T  (ij)  corresponds  to  a  non-interacting  collision  between  two  hard 
sphere  molecules 

Tn (i j )  -  -Z2J dS  I  V±J .3  1 6 3  <?  -S.  )  . 


(2-1 7b) 


v. . *a.  ,<o 
ID  ID 


-r-  -r-  -r  -+•-*- 

Here  v^=v^-Vj#  r^-r^r^  are  the  relative  velocity  and  position  of  the 

two  molecules,  a  is  the  diameter  of  the  molecules  and  (T^  is  the  perihelion 

vector  of  the  collision  between  molecules  i  and  j  (=  vector  from  center  of 

j  to  center  of  i  at  time  of  contact) .  The  rotation  operator  Rn  transforms 

uiD 

■  )■  ■ )  y  ^ 

the  velocities  v, ,  v.  pviov  to  the  collision  into  the  velocities  v! ,  v'. 

ID  i  3 

•j* 

after  the  collision 


+  ■  A 

T  In  this  report  a  symbol  a  always  indicates  the  unit  vector  in  the 
direction  of  the  vector  It. 
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V!  =  R  V.=  V.  —  (V.  .  *8.  .)8.  .  , 

i  a.  .  i  i  13  13  1] 

v!  -  R13  v.=  v.  +  (v.  .*8.  .)8.  .  .  (2“ 

:  3  3  l3  13  13 

The  properties  of  these  binary  collision  operators  were  discussed  in 


(2-18) 


detail  in  previous  reports  {26,27] . 

As  shown  by  McClure  and  Dorfman  [17]  and  also  discussed  by  Kuperman 
[19]  the  time-displacement  operator,  defined  in  (2-12),  can  be  expressed 


in  terms  of  these  binary  collision  operators  as 

N 

T  (1 . N;X)  =  e‘IvIT  11001  . 

t  V— 1 


(2-19) 


where  G„  is  the  same  resolvent  operator  as  used  in  AEDC-TR-72-142  [27] 

N  AT  Ilt 


Gn(1,...N)  = 


h  *  k+Zz!(ij) 


(2-20) 


1=1  1  i<3 

In  the  absence  of  the  object  the  time-displacement  operator  reduces  to 


Tt(l,...N)  =  eN 


(2-21) 


The  explicit  expressions  (2-17)  for  the  binary  collision  operators 
T(ij)  refer  to  a  gas  of  hard  spheres.  However,  the  crucial  assumption 
enabling  us  to  represent  the  time-displacement  operators  by  (2-19)  and 
(2-20)  is  the  assumption  that  the  size  of  the  molecules  is  small  com¬ 
pared  to  the  mean  free  path.  Since  we  shall  apply  the  theory  to  the 
case  that  even  the  macroscopic  object  is  small  compared  to  the  mean 
free  path,  this  assumption  is  satisfied  for  any  gas  of  molecules  with 
finite  interaction  range.  Therefore,  we  may  use  (2-19)  and  (2-20)  to 
represent  the  time-displacement  operator  for  any  gas  of  molecules  in 
the  nearly  free  molecular  flow  regime,  provided  that  the  hard  sphere 
cross  section  in  (2-17)  is  replaced  with  the  actual  cross  section  of 


the  molecules  under  consideration.  For  the  same  reason  we  do  not 
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distinguish  between  T(ij)  and  T (ij)  operators  as  was  needed  in  the  theory 
of  transport  properties  of  moderately  dense  gases  [27] . 

2.3  Cluster  Expansion. 

According  to  (2-4)  and  (2-13)  the  force  E(t)  on  the  object  is 


i(t)  =  -  fr  /  dT  D  ;X;0)  T  (1,...N;X)  Ptf  )  . 

dt  /  N  N  N  t  N 


(2-22) 


Using  the  form  (2-19)  for  the  time-displacement  operator  T  we  thus  obtain 


E(t)=  -  Ar„  Dm(T  ?X;0)  T  [g  +>T(kX)  P 

f  N  N  N  t  N 


(2-23) 


The  operator  G„  is  associated  with  the  streaming  of  the  molecules  without 
any  interactions  with  the  object.  In  the  absence  of  the  object  the  total 
momentum  P  is  conserved/  so  that  G„  and  ?  commute:  G  ?  =  $G„  =  0.  •  Moreover 

N  N  N 


T (kX)p£  =  0  for  £+k  and  (2-23)  becomes 

E(t)  "  -5^^  DN(rN;X70)  Tt(1'  •••»»»x)T(kx)p]c  ,  (2-24) 

which  for  a  gas  of  N  identical  molecules  reduces  to 

E(t)»  -Ny*dTN  DN(r-N;X;0)  Tt(l,...,N;X)T(lX)P1  .  (2-25) 

We  now  write  the  time-displacement  operator  Tt(l,...N;X)  in  the  form 
of  a  cluster  expansion  [17,19]. 


Tt(l,...,N;X)  =  Ut(l)Tt(2,.. 


+ 

+ 


Ut(li)Tt(2, 
.  +ut(i,... 


. . .,  i-l,i+l, 
,N;X) 


,N;X)  +  Ut(l;X)Tt(2, .. .,N)  + 

...N;X)  +  Ut(li;X)Tt(2,...,i-l,i+l, 

(2-26) 


Here  we  are  introducing  Ursell  operators  defined  by 


ut(l)  =Tt(  1)  , 

Ut(l;X)  =  Tt(l;X)-  Tt(l)  , 
ut(12)  =  ^td2)- Tt(l)7*t(2) 

Ufc(12;X)  =  Tt(12?X)-  Tt(l;X)Tt(2)  -  Tt(2;X)Tt(l) 

-  Tt(12)  +  2  T^dJT  (2) 

.  .  .,  etc.  (2-27) 
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We  substitute  the  expansion  C2-26)  for  the  time-displacement  operator 
T  (1,...N;X)  into  the  expression  (2-25)  for  the  drag  force.  Since 
the  T  operators  in  (2-26)  do  not  affect  the  momentum  p^  of  molecule  1, 


we  obtain 


E(t)=  -^1^  DN(rN;x?o)  |(Ut(l)  ♦  Ut(l?X)}  + 

+  (N-l)  (Ut(12)  +Ut(12;X)}+  .  .  .+  1^(1, .  .  .N;X)Jt  (IX)  ,  (2-28) 
where  we  have  again  used  the  fact  that  the  gas  consists  of  N  identical 


molecules .  Introducing  reduced  distribution  functions  defined  as 


Fs(Xi,.,.,Xg;X;t)  -  ^y^Xs^.-.dXjj  DN(rN;Xft)  ,  1 

we  obtain  in  the  thermodynamic  limit  N-**>,  N/T2=n 

E(t)=  -n/dx^  F1(x1;X,-0)  [Ut(l)+Ut (1;X)  IKlX)^  + 

*4  /"  •*» 

-n  1  dX;Ldx2  F2(xlfX2;X?0)  [Ut(12)+  Ut(12fX)J  TUX^  + 


(2-29) 


+  .  .  •  » 

This  expression  relates  the  force  to  the  distribution  functions 


(2-30) 


F  (x  . ..x  ; X; 0)  at  the  initial  state.  Following  McClure  and  Dorfman 
S  X  s 

we  assume  that  the  initial  state  is  unaffected  by  the  presence  of  the 


object 


(2-31) 


F  (x  . . . x  ;X; 0)  =F  (x, . .  .x  ;0)  .  I 

sis  sis 

We  also  introduce  the  usual  assumption  that  the  molecules  are  un¬ 


correlated  initially,  so  that  F  (x, . . .x  ;0)  may  be  approximated  by  a 

sis 

product  of  single-particle  distribution  functions 


(x_  . .  •  x  ,  0)  —  I  F  (x. ;  0)  . 

.  1  s  Ui  1  1 


(2-32) 


The  extent  to  which  this  assumption  is  justified  has  been  discussed  by 
Dorfman  and  Cohen  [28].  In  the  absence  of  the  foreign  object  the  time- 
displacement  operator  (1, , . ,N)  may  be  reversed  in  time;  therefore 
the  single-particle  distribution  function  F^(x^;t)  in  the  absence  of 
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the  object  may  be  written 
F 


(x^t)  =  .  .dx^  T_t  (I# . .  .N)  DnCTn;X;0) 


(2-33) 


or  with  (2-26)  and  (2-29) 


F^x^t)  -  U^U)  FjtXjj 


;  0)  +  n j’  dx2U_ 


t(12)F2(«1,x2;0)  +  .  .(2-34) 


Using  (2-32)  we  can  invert  (2-34)  to  yield  an  expansion  for  F^x^jO) 


F^x.^0)  =  U+t(l)F1(x1;t)  -nU+fc(l) 


.ydx2U_t(12)j1”u+t(i)F1(xi; 


t)  + 

(2-35) 


Substitution  of  (2-31) , (2-32)  and  (2-35)  into  (2-30)  yields  a  density 

■+1 

expansion  for  the  force  E(t) . 

In  this  report  we  consider  the  drag  force  E  in  the  steady  state 

which  corresponds  to  the  limit  t**°.  In  this  limit  the  gas  is  in 

equilibrium  in  the  absence  of  the  object  and  the  single-particle  dis- 

tribution  function  F^x^t)  becomes  independent  of  the  position  and 

time  t  and  may  be  represented  by  a  displaced  Maxwell  distribution 

2 


limF  (x.jt)  =F(v.;V)  =  (2lBikT) 


_  3j  _  m(vj-V) 


1'  i 


2kT 


(2-36) 


In  this  limit  we  obtain  for  the  density  expansion  of  the  force  E  on  the 
object 

TV 

(2-37) 


zX(t)- 


E  =  lim 
t*«  a 

where,  for  the  purposes  of  this  report,  we  consider  only  the  first  two 

41 

terms  E  and  E, . 
o  1 

Eo(t)— ny*dx1  FCv^-V)  {Ut(l)  +  Ut(l?X)  jTdX)^  ,  (2-38a) 

^1tt)=-n2Jdx1dx2  F(vl}V)F(v2;V)  (Ut(12)  +  Ut(12;X)  JTOXJPj  + 
iydx1dx2[Ut(l)U_t(12)F(v1;V)F(v2,-^)]  (Ut(l)  +  Ufc(l;X)  }T(1X) pr(2-38b) 


E 

+n 


t  The  operators  U  (l...s)  do  not  involve  the  object  and  the  corresponding 
time  reversed  operators  U  (l...s)  are  well  defined.  The  operators 
U  (l...s;X)  on  the  other  hand  depend  on  the  stochastic  interaction  of 
the  molecules  with  the  object  and  therefore  cannot  be  reversed  in 
time  [17] . 
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Here  we  use  the  convention  that  operators  inside  the  square  brackets  [  ] 
do  not  operate  on  terms  outside  these  brackets.  In  deriving  (2-38)  use 
is  made  of  the  fact  that  (x)F (v^;V)  =  F(v^;V).  Since  the  adjoint 

operator  (U  (1)U  .(12))^'  =  U  (12)U  (1)  and  since  U  (1)U  (1)  =  1, 

+t  — t  +t  — t  — t  +t 

(2-38b)  reduces  to 

^(t)  =  -n2/dx1dx2  F(v1;V)F(v2;V){Ut(12;X)-Ut(12)U_t(l)Dt(l;X)}T(lX)5i. 

(2-38c) 

2.4  Binary  Collision  Expansion 

->■  "► 

The  various  terms  E  (t)  in  the  expansion  (2-37)  for  the  force  E 

(X 

contain  the  dynamics  of  a+1  molecules  and  the  object.  When  the  density 
of  the  gas  is  sufficiently  small  so  that  the  mean  free  path  is  large 
compared  to  the  size  of  the  molecules  we  may  neglect  the  probability  that 
two  collisions  occur  simultaneously  and  the  dynamical  processes  reduce 
to  sequences  of  successive  collisions  among  the  molecules  and  the  object. 
In  order  to  classify  the  various  collision  sequences  that  contribute  to 
the  drag  it  is  convenient  to  represent  the  time-displacement  operators 
T^(l...s)  and  T^(l...s;X)  in  terms  of  a  binary  collision  expansion. 

A  precise  formulation  of  the  binary  collision  expansion  for  the 
time-displacement  operator  of  a  gas  of  hard  spheres  in  the  absence  of 
an  object  was  developed  by  Ernst  et  al.  [25] .  The  extension  of  this 
procedure  to  the  time-displacement  operator  in  the  presence  of  the  object 
was  formulated  by  McClure  and  Dorfman  [17]  and  further  discussed  by 
Kuperman  [19] .  Following  these  authors  we  note  that  the  time-displacement 
operator  (2-19)  has  the  form 
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T  (1. ..s)  =  e 


t[L°(l...s)+K(l...sjX)] 


where 


L°(l...s)  =  >  r1- 


and 


V  Pi  .  3 

4^  m 

s  s 

K(l...s;X)=  +  ^T  (kX)  . 

i<j  k=l 

The  time-displacement  operator  satisfies  the  relation 

T 


(2-39) 


(2-40) 


(2-41) 


rt  -  S«  +  A  Si-T  K  TX  ' 


(2-42) 


where  S®  is  the  streaming  operator  representing  the  free  streaming  of 
the  molecules  without  any  interactions 

(2-43) 

The  binary  collision  expansion  is  obtained  after  successive  iteration 
of  (2-42) 


« / 1  _ \  _  tL° ( 1 • • • s ) 


S®(l...s)  -  e 


t  t 

=  S°  +  /  dT  S°  K  S°  +  /dT.  / dT  s'  K  S°  K  S' 

t  t  J  1  t-Tl  TX  j  2  t-Tl  Tl-X2  T 


K  S'  + 
2 


(2-44) 


Introducing  the  convolution  product 
t  t 


■/ 


f  *  g  =  IdT  f(T)g(t-T)  = 


= JdT' 


f(t-T')g(T*)  , 


(2-45) 


o  o 

where  f(t)  and  g(t)  are  functions  of  the  time  t,  the  binary  collision 
expansion  (2-44)  may  be  written  as 

T  =  S'  +  S°*KS°  +  S°*KS°*KS'+  .  .  .  (2-46) 

v  t 

The  operator  K  represents  a  collision  either  between  two  molecules 
or  a  molecule  and  the  object  and  the  operator  S®  represents  the  free 
streaming  of  the  molecules  between  collisions.  Thus  the  terms  in  the 
binary  collision  expansion  (2-46)  correspond  to  sequences  of  zero,  one, 
two,  etc.  successive  collisions  among  the  s  molecules  and  the  object. 
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However,  only  those  terms  in  (2-46)  will  be  different  from  zero  which 
correspond  to  sequences  of  collisions  that  are  permitted  by  the  laws 
of  mechanics.  The  number  of  physically  possible  sequences  of  collisions 
depends  on  the  geometry  of  the  object.  In  this  report  we  shall  limit  our 
considerations  to  convex  objects.  A  molecule  emitted  by  a  convex  object 
cannot  return  to  the  surface  of  the  object  unless  it  first  collides  with 
another  molecule,  so  that 

S°*T(iX) S°*T (iX) S°  =  0.  (2-47) 

2 

Furthermore,  the  expansion  parameter  is  the  inverse  Knudsen  number  no  L, 

where  L  measures  the  size  of  the  object.  Terms  in  (2-46)  that  contain 

more  than  one  binary  collision  operator  corresponding  to  a  collision 

between  two  molecules  may  be  neglected  since  they  are  of  higher  order 
3 

in  no  .  A  complete  list  of  the  combinations  of  binary  collision  operators 
that  yield  vanishing  contributions  was  presented  by  Kuperman  [19] . 


Substitution  of  (2-27)  into  (2-38)  yields 

/"  “►  “►  “► 

EQ(t)=  -n/dx-j^  F (v^;V)  T^U-JX^  , 

E1(t)=  -n^*dxidx2  F(vi;V)F(v2;V) 


(2-48a) 


(2-4 8b) 


•  IT  (12;X)-T  (1)T  (2;X) -T  (1;X)T  (2) -T  (12) +2 T  (1 )T  (2) }t (IX) . 

U  t  t  ttt  tt  1 

Introduction  of  the  binary  collision  expansion  (2-46)  into  (2-48)  then 


leads  to  [17,19] 

E  (t) =  -n/dx..  F (v. ;V)T (IX) p. 


with 


1=  -n/dx 


E,(t)=  -n  Idx  dx_  F (v  ;V)F (v  ; V)  (B,+BJ 


3  4'  ' 


(2-49a) 

(2-49b) 


B3=  S°*K(12;X)S0*KC12;X)S°T(1X)P1  ,  (2-50^) 

B4=  S0*K(l-2;X)S0*K(12;X)S0*Ka2;X)S°T(lX)p1  .  (2-50b) 
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For  convex  objects  the  binary  collision  expansion  of  Eq  terminates  after 
the  first  term  at  a  result  of  (2-47)  and  the  binary  collision  expansion 
of  terminates  after  the  second  term  as  shown  by  Kuperman  [191 .  If 
we  substitute  K(12;X)  *  T(12)+T(1X)+T(2X)  into  (2-50)  and  omit  all 
vanishing  termsj  and  reduce  to 

B3=  S0*{t(1X)+T(2X>}S°*T(12)S°T(1X)p1  ,  C2-51a) 

B  *  S°*{T(lX)S0*T(2X)+T(2X)S0*T(lX)}S0*T(12)S°TtlX)p1.  (2-51b) 

Since  we  may  interchange  the  integration  variables  x^  and  x ,  we  adopt 
the  convention  that  molecule  1  is  identified  as  the  molecule  in  the  left 
most  T(iX)  operator,  i.e.  molecule  1  is  the  molecule  that  collides 
initially  with  the  object.  We  may  therefore  replace  (2-51)  with 

B3=  S0*T(lX)S“*T(12)S0{TCLX)4T(2X)Hp1+p2},  (2-52a) 

B4=  S°*T(lX)S0*T(2X)S0*T(12)S°{T(lX)4T(2X)}{p1+p2}.  (2-52b) 

The  term  B3  accounts  for  the  effect  of  sequences  of  three  successive 
collisions  among  two  molecules,  molecules  1  and  2,  and  the  object,  and 
the  term  B^  is  related  to  sequences  of  four  successive  collisions  among 
two  molecules  and  the  object.  A  complete  classification  of  these  collision 
sequences  was  made  by  Kuperman  [19] .  In  order  to  enumerate  these  collision 
sequences,  we  use  (2-16)  and  (2-17)  to  express  B3  and  B^  more  explicitly  as 

B3  =  S°*  T1(1X)S°*T1(12)S°T(1X)P1+  (Rl) 

+  S°*  Tn(lX)S°*T1(12)S°T(lX)p1+  (R2) 

+  S°*  T1(lX)S0*T1(l2)S<,T(2X)p2+  (Cl) 

+  S°*  Tn(lX)S°*T1(12)S°T(2X)p2+  (C2) 

+  S°*  T1(lX)S°*Tn(12)S°T(2X)p2+  (HI) 

+  S°*  Tn(lX)S°*Tn(12)S°T(2X)p2  ,  (H2)  (2^53,) 

and 
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§4  =  S#*T1ClX)S0*T:LC2X)S0*T1a2)S0TClX)p1+  (R1C1) 

+  S°*Tn(lX)S0*T1(2X)S0*T1(12)S°T(lX)p1+  (R2C1) 

+  So*T1ClX)S0*TnC2X)So*T1C12)S0T{lX)p1+  (R1C2) 

+  So*Tn(lX)S°*Tn(2X)So*T1(12)S0T{lXlP1+  (R2C2) 

+  S°*Ti{lX)S°*Ti(2X)S0*T1(12)S0T{2X)p2+  (C1R1) 

+  S°*TnClX)S0*TiC2X)S0*T1C12)S°T(2X)p2+  (C2R1) 

+  So*T1(lX)So*Tn(2X)So*Ti(12)S0T(2X)p2+  (ClR2) 

+  S°*Tn(lX)S0*Tn(2X)S°*T:L{12)S0T(2X)p2+  (C2R2)  .  (2-54) 

According  to  (2-53)  can  be  decomposed  into  a  sum  of  six  terms 
corresponding  to  six  different  types  of  sequences  of  three  successive 
collisions  among  two  molecules  and  the  object.  In  analogy  to  the 
previous  work  of  Sengers  [29,30]  we  refer  to  these  sequences  as  R1  and 
R2  (recollisions) ,  Cl  and  C2  (cyclic  collisions)  and  Hi  and  H2  (hypothe¬ 
tical  collisions) .  These  collision  sequences  are  represented  by  the  six 
diagrams  in  Fig.  2.  In  these  diagrams  the  lines  with  labels  1  and  2  indi¬ 
cate  the  trajectories  of  molecules  1  and  2.  The  vertical  line  represents 
the  object.  The  molecules  traverse  their  trajectories  in  the  direction 
indicated  by  arrows.  An  interacting  collision  between  the  two  molecules 
or  a  molecule  and  the  object  causes  a  change  in  the  direction  of  the 
trajectories.  In  a  non-interacting  collision  the  molecules  continue 
to  proceed  in  the  direction  of  their  original  trajectories;  in  the  diagrams 
we  add  a  shaded  region  wherever  we  want  to  indicate  the  occurrence  of  a 
non-interacting  collision. 

Molecule  1  initially  strikes  the  object  and  is  either  reflected  by 
the  object  (Rl,  Cl,  Hi)  or  passes  through  the  object  (R2,  C2,  H2) .  It 
then  collides  with  molecule  2  such  that  molecule  1  either  collides  with 
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the  object  again  (Rl,  R2) ,  or  it  causes  molecule  2  to  collide  with  the 
object  (Cl,  C2),  or  it  prevents  molecule  2  from  colliding  with  the  object 
(Hi,  H2) .  On  comparing  Pig.  2  with  Fig.  5  of  AEDC-TR-69-68  [30]  we  note 
that  these  collision  sequences  are  closely  analogous  to  the  collision 
sequences  earlier  encountered  in  the  first  density  correction  to  the 
transport  properties  of  a  moderately  dense  gas.  In  fact,  the  new  collision 
sequences  are  obtained  if  in  the  earlier  three-particle  collision  se¬ 
quences  one  of  the  molecules  is  replaced  with  the  object. 

According  to  (2-54)  can  be  decomposed  into  a  sum  of  terms 

corresponding  to  eight  different  types  of  four  successive  collisions 

among  two  molecules  and  the  object.  These  sequences  are  represented 

schematically  by  the  diagrams  of  Fig.  3.  As  noted  earlier,  for  convex 

objects  we  do  not  have  to  consider  any  sequences  of  five  successive 

collisions  among  two  molecules  and  the  object  in  calculating  the  first 

inverse  Knudsen  number  correction  to  the  drag  force. 

•  •  ”►*  *>■ 

The  relative  magnitude  of  B^  and  B^  is  determined  by  the  probability 

that  the  corresponding  collision  sequences  will  occur.  In  evaluating  the 

first  density  correction  to  the  transport  properties  of  a  gas  of  hard 

spheres  Gillespie  and  Sengers  noted  that  the  contribution  from  se- 

-4 

quences  of  four  successive  collisions  is  only  10  times  the  contribution 
from  sequences  of  three  successive  collisions  [31] .  A  preliminary 
numerical  analysis,  conducted  by  Kuperman  indicated  that  for  the  drag 
problem  sequences  of  four  successive  collisions  would  yield  a  correction 
of  less  than  one  percent  [19] .  Therefore,  to  obtain  a  realistic  estimate 
of  the  magnitude  of  the  drag  force  we  shall  only  consider  the  contributions 
from  three  successive  collisions  and  approximate  (2-49b)  by 

E1(t)=  -n 2f  dx1dx2  F(vi;V)F(v2;V)B3  .  (2-55) 
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2.5  Explicit  Formulation  of  Collision  Integrals 

The  drag  force  in  the  free  molecular  flow  regime  is  determined  by 
(2-49a) .  Introducing  the  explicit  form  C2-16)  for  the  binary  collision 
operator  T(1X)  and  replacing  the  momentum  integration  with  a  velocity 
integration,  we  thus  obtain 


E  = 
o 


=  -inn  j  dv^f  (v^V)  J  dA^^  Iv^nJ  j  dv^  n.Cv£  |  V1)  (vj-V^ 


)  r 


C2-56) 


v  A 

v^>n^<o 


with 


f(vi?V)  = 


m 


V ;  - 


HI  tV. -Vi 


27TkT 


2JCT 


C2-571 


The  force  E^,  to  be  added  in  the  nearly  free  molecular  flow  regime 


is  given  by  (2-55)  and  can  be  decomposed  into 

+  +  +  4  4  +  + 

E„  =  +  E _  +  E„,  +  E__  +  E„,  +  E. 


(2-581 


'1  “R1  '  "R2  '  “Cl  '  “C2  '  “HI  ’  “h2  * 

where  each  term  is  uniquely  related  to  the  corresponding  term  in  (2-531 , 

Each  term  in  (2-58)  may  be  written  explicitly  as  a  collision  integral 
in  terms  of  the  initial  velocities,  intermediate  velocities  and  final 
velocities  in  the  collision  sequences  of  Pig.  2  as  shown  by  Kuperman 
{19].  As  an  example  we  consider  the  contribution  ERl.  Using  the 
definition  (2-451  of  the  convolution  productjERl  may  be  written  as 

2  .  f  -4-  -4  4-4- 

lim  J dTjdr  2dv1dv2  f  (v^-V)  f  (v2;V) 


eri  *  -"n 


t-X» 
-T 


*/dr/dT'  S®_T  T1  (IX)  S»_T,  S® ,  TUXlVj.  (2-591 

o  o 

It  is  convenient  to  introduce  the  transformation  T-T'  and 
T2=  T* ,  so  that 


ERr  -”n 


/  4  4  4  4  4  -4  -4  -4 

lim  ydr1dr2dv1dv2  f  (v^-V)  f  (v2;V) 
■|dT2  /dTl  SWVT2)  ^UX)  T1 


(12)  S®  T(lXlyi.  (2-60) 
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The  operators  S° (12)  are  free  streaming  operators,  so  that 


S°r,=  r,  +  v  T 
Til  1 

-*■  -+  -*■ 
S°r  =  r_  +  v.t 
T  2  2  2 


s4,-  \  , 

->■  ->■ 
S°V_“  v„  . 
T  2  2 


(2-61) 


The  left  most  operator  S°  .  .  .  in  (2-60)  translates  the  integrand 

t-tT^  t2) 

according  to  the  free  particle  motion  (2-61) .  Since  the  integration 
extends  over  the  entire  phase  space  of  molecules  1  and  2,  the  integral  is 
invariant  under  such  a  transformation  and  (2-60)  reduces  to 

oo  oo 

-mn2  /dr ^dr 2 dv^dv^ f  (v^  V)  f  (v2 ;  V)  j  ^2f6Xl  {1X)  T±  (12*  ST  T^x^1' 

J  o  o  1  2(2-62) 

In  Fig.  4  we  consider  again  a  typical  recollision.  The  lines 

represent  the  trajectories  of  molecules  1  and  2  and  the  circle  indicates 

the  object.  Molecule  1  collides  with  the  surface  of  the  object  at  the 

•f 

position  (first  collision) ,  it  subsequently  collides  with  molecule  2 

(second  collision)  and  then  it  collides  with  the  object  again  at  the 

position  R2  (third  collision) ,  The  first  collision  may  either  be  an  in- 

teracting  collision  (Fig.  4a)  or  a  non-interacting  collision  (Fig.  4b) . 

We  indicate  the  initial  velocities  of  the  molecules  by  v^,  the  velocities 

after  the  first  collision  by  v' »  v',  the  velocities  after  the  second 

1  * 

collision  by  v£,  v"  and  the  velocities  after  the  third  collision  by  ^ 1 1 , 
v”  1 ;  we  indicate  the  relative  velocities  similarly  by  v^2=  , 


v'  =  v 1  -v ' ,  v' '=  v1 '-v 
12  1  2'  12  1 

recollisions 


2  '  12 


vi  1 1 _vi » i 
1  2 


We  note  that  for  the 


v ' =  v_,  v' 1 '=  v' ' 
2  2  2  2 


and 


(Rl-sequence) 


(2-63a) 


v'=  v, ,  v'=  v„,  v' ' 
1  12  2  2 


(R2-seguence) . 


(2-63b) 


Using  this  notation  and  the  definitions  (2-16),  (2-17)  and  (2-61)  of  the 
operators  T(iX),  T(12)  and  S°(12),  the  expression  for  £  may  be  vuritten 

Rl 
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Figure  4 .  Recollisions 


34 


AEDC-TR-74-79 


CO  00 


t  2  2 

Er1=  -mn  a 


y dr  1dr2dv1dv2  f (v^V)  f  ^jVj^dT.^dT^y' ^1  I  vl'”l I 


6 3 (r-R, ) 


v^»n<o 


/d»l  /d3i2  l7i2'9l2l  6  <?l2+?i2Tl-512> 


Vf  .0  <0 

12  12 


•  |dA2  lv^-^1  VJTJ+  v|'T2-R2)ydv^"n(v|"|v|')  (v”'-v”),  (2-64) 

v  '  '■•n  <o 

X  4 

where  cLA^and  dA^  are  the  surface  elements  at  the  positions  R^  and  R^r 
respectively.  The  time  is  to  be  identified  as  the  time  between  the 
first  and  the  second  collision  and  the  time  T2  as  the  time  between  the 
second  and  the  third  collision.  The  first  two  delta  functions  in  (2-64) 


can  be  readily  integrated  and  we  obtain  [19] 

00  00 

-mnJ“ 


£..=  2 


R1 


02Jdv1dv2  f  (v^V)  f  <V2iV)J dT^|dT„ 


K  IVSil  /d?l  n(’i|,'i)/a8i2l?i,2'ai2l /“a  l?r-8: 


'r¥° 


v  . 1  .  <o 

12  °12 


v» 

1  2 


-a  ->■  -*■ 

,<s  <VV 


j+  v^T1+  v^’T2) j dv”'n(vj"|v|')  (v^'-v”).  (2-65a) 

Taking  cognizance  of  the  difference  between  T1(1X)  and  Tn(lX),  we  obtain 


for  E 


R2 


00  00 


2  2/  ■>  -*  -*  ->  /  r 

=  +mn  tr/dv1dv2  f  (v^VJf  (v^Vj/dTj* 

*K  I^L*“ll  /dSl  2  I  ^12  ’  612 1  f**2  l^'‘*“21  {3(VVVl+;i 


V 


v^*n^<o 


v  •$  <o 
12  12 


v-.572<o 


•Jfe* '  ' Tl  (v| '  1 1  1  ' 4j ' )  .  (2-65b) 

Following  the  same  procedure  we  obtain  for  the  contributions  EC1 


and  from  cyclic  collisions 
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oo  <x> 


SQ=  -mn2C^yrdv1dv2  f  (v^)  f  (v^V)  j&T^dz 

’/j*l  |v1-n1|/dv'n(v*|vi)/d012  |v'2-®12|  fd 
^  •h«<o  v'.  <o  ±i  »„S 


dA  v''*n_ 
2  1  2  21 


1  1 


12  12 


^2  V° 


•63(R1-R2+  vjr^  v^*T2)Jdv^,,n(v*,,|v”)  (v^'-v^') 


t2-66a> 


E  =  +mn2<72/dv  dv  f(v, ;V)f(v 


y 


1  2 


OO  00 

fdT. 


jV)  JdTl/d 

o  o 

/*!  I  V^l I  f*  12  l’l2-ai2l/“2  ‘,'VV;iY;i'V 


v1*n1<o 


v  <o 

12  12 


v' ' *fi  <o 
2  2 


,  j dv2 ' 1  f|  ( '  '  |v^')  (v^'  '-v*') 


where  the  various  velocities  are  now  indicated  in  Pig.  5  for  which 


•  •  •  =  »'■ 


and 


v  =  v  «  V 
2  2  1 


v,,  v'=  v  .  v 


i"_  „n 


(.Cl-sequence)  , 


(C2-sequence)  . 


{2-66b) 


(2-67a) 


(2- 67b) 


The  contributions  E  _ ,  and  E  _  from  hypothetical  collisions  are 

Hi  HZ 

00  00 

EH1=  -Hnn2CT2y’dv1dv2  f  (v^-V)  f  (v2;V)^  dT1^dT2 

’  ydvjrUvjJv^  yd812|v^2-912[  |dX2  tv2*Ji2| 


v^*ft^<o 


Ti2-812<0 


V  *ft_<0 
2  2 


•  63  (R1-R2+  vJTj  +  V2T2) y  dv^'  ’U(v”  '  |v”  )  (^'  '-v^') 

do  TO 

2  2/+  +  f  [ 

2=  -nm  a  I d»1dv2  f  (v^-v)  f  (v2>-v) fd^  /  dt2 

/**!  1*1*1  1^12  I^12*ei2l  l^2'a2l  63(WVl+  V2) 


t2-68a) 


vj^-ft^o  ^12  "^12<0 


v  -n  <o 
2  2 


*y dv”  »T1  •  | v^ • )  (v*  " -v  ” ) 


(2-68b) 
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where  the  various  velocities  are  now  indicated  in  Fig.  6,  for  which 

V2,=V2=V2'  Vi',=V3.,=Vi  (Hl-sequence)  ,  (2-69a) 

and 

v''=v'=v  ,  v' 1  '=v'  '=v'=v  CH2-sequencel  .  (2’~69bi 

/  ^  ^  i.  X  1  Ji 

If  the  molecules  of  the  gas  are  approximated  by  hard  spheres  a 
major  simplification  occurs,  since  we  can  prove  that  the  recollisions 
and  the  cyclic  collisions  yield  identical  contributions 


Rl 


■y  -y 
E  =E 
C2  R2 


C2-701 


This  theorem  is  proved  in  the  Appendix. 

In  conclusion,  we  may  approximate  the  first  inverse  Knudsen  number 
correction  to  the  drag  force  by 

V^l  +  %2  +  2<JRl+SR2)  •  C2‘71) 

with  E  and  E  given  by  (2-65)  and  E  .  and  E  given  by  C2-68) . 

Rl  R2  HI  Hz 

The  term  E  represents  a  "loss"  term  which  accounts  for  the  fact  that 

Ill 

molecules  reflected  from  the  object  prevent  some  molecules  of  the  incident 

beam  from  striking  the  object.  The  term  2E^  represents  a  "gain”  term 

which  accounts  for  the  fact  that  molecules  reflected  from  the  object 

cause  some  additional  molecules  to  collide  with  the  object.  The  terms 

Eh2  and  2Er2  account  for  a  perturbation  of  the  velocity  distribution  as 

a  result  of  the  fact  that  the  presence  of  the  object  leads  to  a  region 

that  is  inaccessible  to  some  gas  molecules.  The  terms  E  and  2E  _ 

H2  R2 

vanish  in  the  infinite  Mach  number  limit  [19] ,  but  need  to  be  included 
in  determining  the  drag  force  at  small  speed  ratios. 
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CHAPTER  III 


DRAG  COEFFICIENT  OF  A  DISC  IN  THE  NEARLY  FREE 
MOLECULAR  FLOW  REGIME 


3 . 1  Introduction 


As  a  first  application  of  our  method  we  consider  the  drag  force 
exerted  on  a  disc.  A  study  of  the  drag  of  a  disc  became  very  appealing, 
since,  just  prior  to  this  research  effort,  Willis  had  obtained  an 
analytic  solution  for  this  case  in  the  infinite  .Mach  number  limit  [20,21]  + 
The  situation  is  illustrated  in  Fig,  7.  A  disc  with  radius  R  is 
placed  perpendicular  to  a  gas  stream  with  flow  velocity  The  disc  is 
located  in  the  XY-plane  and  the  Z-axis  is  taken  in  the  direction  of  -V. 

The  distribution  of  the  velocities  of  the  molecules  in  the  gas 
stream  can  be  represented  by  a  displaced  Maxwell  distribution  (2-57) 


(3-1) 


The  molecules  are  assumed  to  be  reflected  diffusively  at  the  surface 


so  that  the  transition  probability  n(v'  fv)  in  (.2-15)  is  given  by 

mv|Z 

Tl(v'|v)  =  n(v')  =2~  (v'-n)e  2kT  0(v’-n)  ,  (3-2) 


+  We  are  indebted  to  Dr.  A.  G.  Keel  of  the  Naval  Ordnance  Laboratory 
for  informing  us  about  this  work  prior  to  publication. 
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where  0  is  a  Heaviside  function  defined  as 

0(x)=  1  for  x  >  o  ,  (3-3) 

0  (x)  =  0  for  x  <  o  . 

The  velocities  v  and  v'  are  the  velocities  of  a  molecule  before  and  after 
striking  the  object.  The  vector  n  is  the  normal  vector  of  the  object 
surface.  At  the  upper  side  of  the  disc  surface  n  is  the  unit  vector  in 
the  +Z-direction;  at  the  lower  side  of  the  disc  surface  ft  is  the  unit 
vector  in  the  -Z -direction.  The  temperature  T  in  (3-2)  is  the  temp¬ 
erature  of  the  object  which  is  taken  to  be  equal  to  the  temperature  T 
of  the  gas  in  (.3-1)  . 

In  order  to  exhibit  the  dependence  of  the  drag  force  on  the  inverse 

Knudsen  number  explicitly  it  is  convenient  to  introduce  dimensionless 

quantities.  For  this  purpose  we  measure  all  distances  in  terms  of  the 

radius  R  of  the  disc  and  all  velocities  in  terms  of  the  thermal  velocity 
1/2 

(2kT/m)  of  the  molecules.  We  thus  define 

i  i\2kT/  •  ®  v^2kTj  '  R\m  j  .  (3-4) 

The  magnitude  of  the  dimensionless  flow  velocity  ll  is  the  speed  ratio 

earlier  introduced  in  (1-5) .  The  surface  element  dit  of  the  disc  surface 

can  be  written  as 
2  -*■ 

dA=  R  dr  ,  (3-5) 

where  r  is  a  two  dimensional  vector  in  the  plane  of  the  disc  such  that 
o<r<l.  In  addition  we  introduce  a  dimensionless  distribution  function 


f*(w^;S)  and  a  dimensionless  transition  probability  H * (w!)  defined  as 


f*(w.;S)= 

1  TT/z 


e“  (w, -S) 


C3-61 


42 


AEDC-TR -74-79 


n*(wp=  ^  (w|‘n)e  wi  0(w|‘n)  , 


(3-7) 


/-*-  ->■  r  ■+  -> 

dw^  f*  (w^;S)=  1  and  I  dw|  n* (w|)«  1. 

*  l  2  2 

The  kinetic  energy  of  the  gas  stream  is  UR“  •jnm.V  TO  ,  so 
that/  in  accordance  with  (1-1) /  the  drag  coefficient  of  a  disc  with 
radius  R  is  defined  as 


2E 

D  nrov2TO2  •  (3-8) 

As  a  result  of  symmetry  the  drag  force  has  only  a  non-vanishing  component 

in  the  direction  of  the  flow  velocity  so  that 
-►*  A 

E=  E-V  .  (3-9) 


In  this  report  we  define  the  Knudsen  number  K  as  the  ratio  of  the  mean 
free  path  ^  over  the  radius  R  of  the  disc.  Since  for  a  gas  of  hard 
spheres  ^  [3] ,  we  thus  define 


K  1=  vT*  irn02R. 


(3-10) 


The  drag  coefficient  CD  in  the  nearly  free  molecular  flow  regime  can  be 
written  as 


C  =  C  +  C.K_1+.  .  .  .  (3-11) 

u  o  i 

The  drag  coefficient  C  in  the  free  molecular  flow  limit  K  1-*o  follows 

o 

from  (2-56) 

Co=~Jp  f*(w;S)Jdr  |w-  fi|  Jdw'T)* (w' )  (w'-W)  .fe.  (3-12) 

w»ri<o 


The  coefficient  of  the  first  inverse  Knudsen  number  correction  can 
be  decomposed  in  analogy  to  (2-71) 

V  CH1+  CH2+  CR2>  '  «-13> 
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Table  I. 


Collision  integrals  for  the  first  inverse  Knudsen  number  correction  to  the 
drag  coefficient  of  a  disc. 


-►  -  _(w  -s)z 

Note:  f*(w^«S)=  TT  e  i  , 

T]*(w|)=  wi  © •  SI)  (diffusive  reflection). 

For  a  sphere:  n1=  r2»  n2« 
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The  explicit  expressions  for  the  collision  integrals  C  . ,  C  ,  C_.. ,  and 

HI  H2  R1 

follow  from  (2-68)  and  (2-65)  and  are  listed  in  Table  I.  The  velocities 
in  the  integrands  refer  again  to  the  hypothetical  collisions  and  recol¬ 
lisions  which,  for  a  disc,  are  shown  schematically  in  Figs.  8  and  9. 

We  note  that  all  collision  integrals  contain  a  factor  of  the  type 
r  -►  ■> 

dw,Ti*(wl )  (w*-w)  which  for  diffusive  reflection  reduces  to 

-w'*  **■ 


/■ 

Idw'n*  (w1 )  (w'-w)=  ld$' e  (w'-n)  (w'-w)=  n-  w  . 


(3-14) 


w'-n>0 


The  drag  coefficient  CD  is  a  function  of  the  speed  ratio  S.  For  a 
disc  it  may  be  decomposed  into  a  contribution  accounting  for  the  force 

exerted  at  the  upper  side  of  the  disc  surface  and  a  contribution 
accounting  for  the  force  exerted  at  the  lower  side  of  the  disc  surface 


c  =  c  ■*■+  C  -  . 

D  D  D 

The  two  contributions  are  interrelated  by 


(3-15) 


C~  (+S)=  -C*(-S) 


(3-16) 


3.2  Drag  coefficient  C^of  a  Disc  in  the  Free  Molecular  Flow  Regime. 

The  drag  coefficient  C  of  the  disc  in  the  free  molecular  flow 
3  o 

regime  is  given  by  (3-12)  which  can  be  integrated  over  the  area  of  the 


disc.  The  drag  coefficient  CQ  of  a  disc  is  thus  the  same  as  the  co¬ 
efficient  Cq  of  a  flat  plate  of  any  shape.  We  decompose  Cq  in  analogy  to 


(3-15)  and  (3-16) 

c  05)  =  C  +(S)+  c  "(S)=  C  +(S)-C  +(-S)  , 
o  o  o  o  o 

and  note  that 


(3-17) 
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Figure  8.  Hypothetical  collisions  for  a  disc. 
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(o)  Rl- SEQUENCE 


(b)  R2- SEQUENCE 

Figure  9.  Recollisions  for  a  disc 
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(3-18) 


4.,^.  2  r+  2 . 2  'fw 

3  «>■  Jd»  «  (V  *  V  - 

w  <o 

f  2  2  2 

|Se  S  +  *^(1+2S2)  (1+erf  S)  +  ^  {e  S  +Sv*tr(l+erf  S) }  |  , 


VTfS- 


where  erf  S  is  the  error  function  defined  as 


erf  s 


(3-19) 


We  thus  obtain  for  the  total  free  molecular  drag  coefficient  C 


2 

C  (S)*—  fe"S  +  m~  +  S)erf  S  +  5 

°  sm  2S  2 


(3-20) 


This  expression  is  in  agreement  with  the  result  obtained  by  previous 
investigators  [32] ,  but  it  has  not  been  reproduced  correctly  in  some 
more  popular  reviews  [1,2/10]. 

In  the  low  and  high  speed  limits  the  drag  coefficient  C  (S)  reduces 

o 

to 


and 


lim  C  (S)  = 
9+o  ° 


(it 1\  I 
fn  s 


4.029 

S 


(3-21) 


lim  C  (S)=  2  .  -  (3-22) 

&*» 

Thus  at  low  speed  ratios  the  drag  coefficient  Cq  (S)  becomes  inversely 
proportional  to  the  speed  ratio  S,  while  in  the  high  speed  limit  the 
drag  coefficient  Cq(S)  becomes  independent  of  the  speed  ratio  S.  On 
comparing  with  (3-8)  we  note  that  this  result  means  that  at  low  velocities 
the  force  is  proportional  to  the  stream  velocity,  as  is  to  be  expected 
[16],  while  at  high  velocities  the  force  increases  with  the  square  of  the 
stream  velocity.  r- 
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3.3  Reduction  of  Collision  Integrals  for 

The  coefficient  of  the  first  inverse  Knudsen  number  correction 
to  the  drag  coefficient  is  determined  by  the  collision  integrals  in  Table  I. 
Since 

cl(s)=  C1+{S)- ci+(-s),  (3-23) 

in  accordance  with  (3-15)  and  (3-16)  ,  we  only  need  to  analyze  the  collision 
integrals  for 

V  C„r<^  +  2fCa  *  CR2>-  (3-24) 

These  collision  integrals  correspond  to  the  case  that  the  normal  vector 
A2  is  in  the  +  Z-direction,  i.e.  in  the  direction  of  -S. 


It  is  convenient  to  introduce  a  set  of  auxiliary  vectors  defined  as 


-*• 

w  , 

= 

W1  - 

W-  w'  w„ 

_ 2_  =  _1  -  2 

(3-25a) 

HI 

■4-  A 

Wi*ni 

w2-ft2  w 

*  w 

lz  2z 

->■ 

w„„ 

— 

-*■ 

W1 

■4 

-  2  = 

■4*  -4- 

w_  w 

f 

(3-25b) 

H2 

l*i*ll 

-»■  . 
w  •  n 

2  2 

w_  w. 

lz  2z 

"hi 

_ 

w1 

"l 

w'  ' 

-Zi _ 

w'  w'  • 

=  1  -  W1  , 

(3-25c) 

w* 

1  1 

w' '  •  n 

1  2 

w'  w'  * 

lz  lz 

-4- 

W1 

w'  ' 
_ 

+1  1 

W1  W1 

WR2 

'  Wi"a2 

=  w  “  wV  • 
lz  lz 

( 3-25d) 

These  auxiliary  vectors  are  two-dimensional  vectors  located  in  the  plane 
of  the  disc.  The  normal  vector  is  directed  in  the  positive  Z-direction 
for  the  integrals  associated  with  the  HI-  and  Rl-  sequences  and  in  the 
negative  Z-direction  for  the  integrals  associated  with  the  H2-  and  R2- 
sequences.  For  future  reference  we  also  define 
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) 


) 


) 


) 


w 

HI 

rm^il3 

A 

WH1  ri  , 

(3-26a) 

FT")2 
'hi  1 

whir5hix  V 

WH1 

WH2 

5H2-fl  , 

(3- 26b) 

(®H2x 

wa2«H2 *  *i>2 

WH2 

WR1 

I^Rl'  *l|! 

5R1‘*1 

(3-26c) 

(5Rix  r/ 

WRlSu  X  V 

“Si  ■ 

to 

“R2-  fl  . 

(3-26d) 

<«E2  x  V2 

“S2«R2  *  81>2 

“L 

As  an  example  we  consider  the  collision  integral  C  *  associated 

£1^ 

with  the  Hl-sequences.  The  6-function  in  (1-1)  may  be  written 

IV"2l63(VV"iTl4V2)  =  6|T2+  •  C3"27) 


Because  of  symmetry  we  can  integrate  over  the  azimuthal  angle  of  the 


position  vector  r^  of  the  first  collision  and  take  r^  in  a  fixed  direction. 


say  in  the  positive  X-direction  as  indicated  in  Figs.  8  and  9.  Using 
id  in 

2V2 


(3-14)  and  integrating  over  T*  and  r^  we  thus  obtain 


'Hi 


(S)  =  — 


r  j  dw^dw^f*  (w1;'S)  f *  (w2;^)  J  Jdr iri  lwlzl 

'Jdh  n*%)  Jddl2l^i2'ei2l  "W2z)  G(_Wlz)  0(_W2z 


C3-28) 


W  *  *Q  <0 
12  12 


with  the  auxiliary  condition 

r  —  Ir  +  (w'*n  )t*w  |<  1  , 

2  1  1  1  1  1H11" 

This  auxiliary  condition  implies 


-W 


u2  w1 
HI  1  1 


C3-29) 


(3-30) 
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The  time  T*  is  to  be  interpreted  as  the  maximum  possible  time  interval 
HI 

between  the  first  and  the  second  collision  that  leads  to  an  Hl-sequence  for 


given  and  w^.  It  then  follows  that 

l 


A;/driri  0(THi"Tr"  -= 


3w'  *h 
1  1 


I(WH1>' 


(3-31) 


with  I (w )  defined  by  (3-26a) .  Integration  over  d  yields 
HI  12 

fad 12  l"lV8l2l  e<““i2-912)  ■  "'ll  ' 

and  the  collision  integral  (3-28)  reduces  to 


(3-32) 


:Hl(S)=  "  3^s2 J i^id”2  f*<w1;S)f*(w2;S)  J dw^e~wl  |wlz| 


(3-33a) 


VF 


'  "2Z)I<WH1>  0(‘Wlz)  e<-”2Z)  0<"L>- 
The  other  collision  integrals  can  be  treated  in  the  same  manner  with 

the  result 

f*  "2ZmiH2)  0tV  0(-W2z)' 

(3-3 3b) 


CH2(S>-  +  if 


f  *  ^2 '  ^  W 


C  "*■  (S)  ■  +  ^ 
RllS)  3TTS4 


fa 


dw2  f* (wi;f) f* (w2;S) 


J dw  je  |  w 


¥  I d»12  <vf  -  0<-»1.>  0<»:.> 0 0(-«, 


lz' 


R1 


lz' 


lz‘ 

* 

lz' 


lz'  12  12 

(3-33c) 


C«(S>=  -  #|  dV»2  £* 

,c 


J dw^dw2  f*  (w^-S)  f*  (”2;^)~^*d^12  I” 


12  d12^ 


<  V  "iz^^  0(Wlz}  0(-Wlz)  0(-W12*ai2)* 


(3-33d) 

The  total  contribution  of  the  four  collision  sequences  to  the  co¬ 


efficient  C1  is  then  determined  by 
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cm(s)*  ch1(s1'  cnlts!  ' 

(3-34a) 

CH2(S>-  CH2lS'-  CH2(-S>  ' 

(3-34b) 

CR1(S>-  CR1(S)-  St''51  ’ 

(3-34c) 

C  (S)=  C  +(S)-  C  +(-S)  . 

R2^  '  R2V 1  R2V  1 

(3-34d) 

a  computational  point  of  view  the  major 

difference  between  C  and 

HJ. 

C„  „  on  the  one  hand  and  C  ,  and  C  „  on  the  other  hand  is  that  in  the 
H2  Rl  R2 

recollision  integrals  the  integration  over  S^2  is  to  be  done  numerically, 
since  the  integrand  depends  on  6^2  via 


wj- «j- 

tin  C^) 

(3-35a) 

Wl=  "l"  <W12’®12,912‘ 

(in  CR2) 

(3-35b) 

3.4  Drag  Coefficient  of  a  Disc  at  High  Speed  Ratios. 

The  calculation  of  the  collision  integrals  for  the  coefficient 
becomes  particularly  simple  when  the  disc  is  placed  in  a  beam  of  molecules 


that  move  with  uniform  velocity  V.  This  limit  is  approximated  when  the 

velocity  V  of  the  gas  stream  is  large  compared  to  the  thermal  velocity 
\J 

(2kT/m)7a  ,  i.e.  S  »  1.  We  may  refer  to  this  limit  as  the  high  speed  limit 
or  Ccold  wall)  beam  limit. 

In  the  beam  limit  all  incident  molecules  are  moving  parallel,  so 

that  two  molecules  will  never  collide  with  each  other  unless  at  least 

one  of  them  first  interacts  with  the  object.  The  contributions  jjim  cH2tS) 

and  j.im  CR2CS),  therefore,  vanish  in  the  beam  limit.  Moreover,  in  the 

beam  limit  the  drag  force  completely  originates  from  collisions  at  the 

upper  surface  of  the  disc,  so  that 

lim.C. (S)=  lim  C„t(S)+  2  lim  C +  (S) .  C3-36) 

&4C0  1  H1  S'**  Rl 
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The  collision  integrals  for  the  drag  coefficient  in  the  beam  limit 

->■ 

are  obtained  when  the  distribution  function  f * (w^;S)  is  approximated  by 
a  6-function 


f*(wi;S)=  6 3 (w^-S) . 


(3-37) 

Noting  that  w^^,  w2=  w^2~  ^S,  wl^=  (S’^12^12  we  tllus  obta:!-n  froiB 

(3-33a)  and  (3-33c) 


as  cHi(s)='ir  sfA le"Wl2  u"hi>  ' 


with 


->i 

W1  -V 

w  =  —*■  +  la  / 

1  Til  " 


HI  w 


(3-38a) 


C3-38b) 


lz 


and 


^4scri(s,'  +  #s 


iJaS’e-"  i2J 


0CW1Z>  S(M12>  ’ 
C3-39a) 


with 


-»-i 


(3-39b) 


These  integrals  can  be  readily  evaluated  analytically.  For  this 

purpose  we  introduce  the  coordinate  transformation 

dw'=  w'4dw'  wdw.,d<{)  , 

1  lz  lz  V  V 

where  w  =  w  in  (3-38)  and  w  =  w  .  in  (3-39)  and  whera  $  is  the  azimuthal 

V  HI  V  RX 

angle  of  the  two  dimensional  vector  w^  with  the  X-axis  as  initial  axis. 

In  terms  of  these  variables  the  quantities  I(w  )  and  I(w  , ),  defined  in 

HI  RX 

(3-26) ,  reduce  to 

The  perihelion  vector  0^^  may  be  expressed  in  terms  of  its  polar  angle 
By  with  the  Z-axis  and  its  azimuthal  angle  4>g  with  the  XZ-plane  as  in¬ 
itial  plane.  The  collision  integrals  then  become 
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e-wU<wHl+1> 


lim  C  (S)  =  -  z^=r-  S  f dw'  w'2  fdw 

g)m  Hi  3TT  J  lz  hi  Hi 

0  0 
2TT  3  , 

•  /<*  -  “*»]  ■  -  3 vf s- 


(3-41a) 


4/J  r0  2  r*  r  P 

CR1(S)=  +  3^  s/dwizwiz/dwsiJd^Jd8a  sin8alcos39al/d<f>1 


Tp4 


“4'Wiz<  4l  +  2wR1tg9aoo3«t-4>0)+  tg20o+3}J  “  <=os*]  =  £§/ 


(3-41b) 


The  total  value  of  the  coefficient  is  then  obtained  from  (3-36) .  We 
thus  conclude  that  in  the  limit  of  high  speed  ratios 


(3-42a) 

(3-42b) 

(3-42c) 


in  agreement  with  the  result  earlier  obtained  by  Willis  [20,21] . 

3.5  Drag  Coefficient  of  a  Disc  at  Low  Speed  Ratios 

At  low  speed  ratios  we  can  expand  the  Maxwell  distribution  (3-6)  in 
terms  of  a  Taylor  series  around  S=0.  If  we  neglect  terms  of  higher  order 
in  S ,  then 

f*(wlfS)f*(w2;S)=  ±3e  VW1  W2  J{l-2S(wlz+  w^)  >  ,  (3-43) 

to  be  substituted  into  the  collision  integrals  (3-33) .  The  zeroth  order 
term  in  (3-43)  does  not  contribute  to  the  drag  coefficient  as  a  result 
of  (3-34) .  We  also  note  that  we  can  integrate  over  w^  in  the  collision 
integrals  for  C  .  and  C_. .  In  the  limit  of  low  velocity  we  thus  obtain 
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from  (3-33)  and  (3-34) 

8  V2  f  -w2 


iim  CH1(S)-  “3lr3S 
&>0 


8  /2" 

y;  CH2(S)-  +iFs 

s-*-o 


j ^  J ^e-wi2wl2(^  -  w2zf  I  twH1)  0(-w2z>  6(W' 

j  d^e""! J, 


lz 
(3-44a) 


)  , 


dw2e“w2(-Wlz-w2z)w12(l^-  w2a)I(wH2) 


■0(wlz)  0(-w2z)  , 


(3-44b) 


lira  C  (S)  =  +3pF 
S-K)  R1  3TT  S 


iVj 


{t-  W2z)  (t*  wii)I‘5Bl - lz'  ”12  -12 - lz- 


,)  0(w,'  )  0(-w'  *9  )  0 (-w’  ' ) ,  (3-44c) 


!!*>  CR2<S)=  -3ft| 
S*0 


Ifarifa-i  f/a»12  I512‘912I 

-*lz-2I>  (4  -i;)  I(V  ®("lz>  0(-*12-»12>  ec— "I.C3-44SJ 


We  recall  that  we  work.  In  a,  coordinate  system  with  the  z-axis  in  the 
direction  of  -S  and  the  X-axis  in  the  plane  through  S  and  r^.  The  polar 
and  arimuthal  angles  of  any  vector  a  in  this  coordinate  system  are  indi- 

A 

cated  by  0  and  <|>  .  We  express  the  perihelion  vector  cr  in  terms  of  a 

3  £1 

polar  angle  0^  and  an  arimuthal  angle  <J>^  in  an  auxiliary  coordinate  system 

•>*  4 

X'Y’Z'  with  the  Z'-axis  in  the  direction  of  w^  (.“w^  ^3-44d) )  and  with 

4-  4- 

the  X ' —axis  in  the  plane  through  S  and  w|2<  The  Cartesian  components  of 

4- 

any  vector  a  in  the  auxiliary  system  are  related  to  its  components  in  the 
original  XYZ  system  by 


/  a  A 

/a 

*  1 

x 

V 

=  R  (0  ,  „,)R  (<j)  ,„,) 
y  wl2  z  Ywl2* 

a 

y 

\az«i 

l  v 

(3-45) 


Here  9w12i  an<*  ^wi2'  are  P°-*-ar  coordinates  of  w|2  in  the  XYZ  system 
and  R  and  R  are  the  same  rotation  matrices  as  those  used  in  earlier 

y  z 

reports  [31,33]. 
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We  conclude  that  the  drag  coefficient  in  the  limit  of  low  velocity 


can  be  written  as 


lim  c  (S)  = 
S+o  1 


lia  C  <S>  +  C  (S)  ♦  2{c  (S)  ♦  C  <S)> 

S+O 


(3-46) 


with 


a/2 

IS  CH1(S)  "  “  3F? 

o  >0 


TT 


2 

p 

„  ,2 

* 

dw  w  e  w2 

j  1  12  -WJ 
dw' w '  e  1 

d8  _sin0  „ 

J 

1  1 

*tt/  ' 

w2  w2 

. 

w 


2TT 

d<j> 


w2 


d0wl,slnewl'  a4„l 

I 

o 


•“h(4  -  W2z)  1  <»H1,  ' 


IT 


a/2 

CH2ts)  =  +  3ir^S 

&>o 


.  2 

* 

.  2 

► 

dw  w  e_W2 

dw  w.  e 

d0  _sin0  _ 

2  2 

1  1 

w2  w2 

• 

2ir 

d<J> 


w2 


d0  _sin0  , 
wl  wl 


2TT 


a*.i< 


% 

/F 


.  (/a  \ 

-w,  -w„  )w,  J—r  -  w_  J 
lz  2z  12^  2  2z/ 


2  "  W2zJI(WH2)  ' 


(3-46a) 


(3-46b) 


pTo 

a“CRllS)  -  +  5f5? 

S-»X) 


00 

0  2 

00 

.2 

TT 

1 

dw  w  e  W2 

dw' w'  e  W1 

d0  _sin0  _ 

2  2 

1  1 

w2  w2 

0 

J 

O 

4 

0 

21T 

d«J) 


w2 


a0wl-5ln6wl' 


2TT 


dd>  — 
VWl'  TT 


TT 


O' 


d0'sin0'  cos0' I 
a  cr  o' 


,2* 

d*c 


% 


«i2  (4  -  ”22)  (4  -  ”ii) 


(3-4 6c) 


CR2(S)  =  '  3?f 

S"*0 


00 

00 

TT 

2 

2 

► 

,  2  -w9 

dw_w_  e  * 

,  2  -w*. 

dWlWl  e  ± 

d8  _sin0  _ 

2  2 

w2  w2 

< 

• 

i 

2TT 

d(f) 


w2 


■% 

d0  ,sin8  . 
wl  wl 


2TT  .TT 


a*„if 


d0 * sinGi I cos0 ' I 
(7  0 '  o' 


2TT 


d*i<-"irW2z)W12 


% 


(4 


R2,B(-ii> 

(3-46d) 
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TABLE  II 


Collision  integrals  for  the  drag  coefficient  of  a  disc  in  the  low 
velocity  limit. 


These  collision  integrals  were  evaluated  numerically  by  a  Monte  Carlo 
method.  That  is,  the  integrals  were  estimated  by  ayeraging  the  integrand 
over  a  set  of  N  random  points  selected  in  the  integration  region  according 
to  a  suitable  predetermined  probability  density  function  J34J ,  For  this 
purpose  we  employed  the  same  method  and  subroutines  that  were  previously 
used  by  Gillespie  and  Sengers  in  a  calculation  of  three-particle  collision 
integrals  for  a  gas  of  hard  spheres  and  reported  in  AEDC-TR-73-171  J31J , 
The  numerical  results,  together  with  their  estimated  standard  deviations, 
obtained  by  averaging  the  integrand  over  50,000  points,  are  presented  in 
Table  II. 
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3.6  Drag  Coefficient  of  a  Disc  at  Arbitrary  Speed  Ratios 

For  arbitrary  values  of  the  speed  ratio  S  we  need  to  substitute  the 
complete  Maxwell  distribution  (3-6)  Into  the  collision  integrals  (3-33) . 
Thus 

Sa151'  -M*  J(sl  / dv"<W2'J)  7aV"i2  »h(f  -  w2Z>i(5hi> 


•0(-w22)  eiw’z>  , 


(3-47a) 


^  2JT_  r+  -(w -s)2  -cw,-s)2__  /r 


C„2<S>-  +3F?/a»1e"'”1'S>7 dV  “12  (T  -  »2z>1(“„2> 

•0(wlz)  0(-w2z)  , 


(3-47b) 


ms) y <£2.-M>2y dt,le-»;2  i ja&i2  u-ej 


(3-47c) 


„  +/0,  2/2 
CR2(S>”  -3F? 

IfH 


ia  i;12-e12f 


T"  WU)I(WP5)  0(W1Z>  0(-Wl5-ai7>  r 


where 


lz'  v  R2 
--*■  *►.  2 


12  12' 
2 


lz' 


■r  t  2  2  V 

J(S)=  ~  JdWj^e  lwlzl  ©(“w^)*  e  S  +  ^  2S(l+erf  S) . 


(3-47d) 


(3-48) 


We  can  express  these  integrals  in  terms  of  the  same  variables  intro¬ 


duced  in  the  preceding  section  and  obtain 

HI 


c"’(s)  =  -  WJ(S1 


dw2W22e’W2 


I 


fTT 


2  -w’ 2 

dw^w^  e  W1  |d8  „sin8 


w2  w2 


.2  TT 


d<f> 


w2 


% 

d8wl-sin6wl' 


!K 


27T 

-S(S+2w_  )  ,  /Vrr  \  ■+ 

dVe  22  wi2  (—  -  w2z/  1  <”„!>"  »-«a) 
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C +  CS)  =  +  -2-^- 
«2  31T3S2 


00 

_  2 

.00 

_  2 

rTT  , 

2  —Wo 
dw„w„  e  * 

dw  w.  e  W^ 

d0  .sin0  . 

2  2 

. 

1  1 

■ 

w2  w2 

r\ 

2T7 

d$> 


w2 


0  0 
%  .2TT 


d0  , sin8  _ 
wl  wl 


-2S(S+w,  +w.  ) 
dd>  ,e  lz  2z 

wl 


*12  -  ”2 J 


Kwh2)  . 


(3-49b) 


-00 


<i(s>  - +  w  J(s> 


dw.w  2e  W2 


d0  sin0  , , 
wl 1  wl 1 


0 

27T  T 


^Wl’lT 


0eo 

2 

2  —  w 1 
dw'w^  e  1 

d8  „sin0  _ 

4 

w2  w2 

< 

0 

2TT 

0 

2ir 


d8^sin0^| cos 0^| 


% 


...  -S (S+w_  )  , 
d(j>^e  2z  wj^ 


«-ii) 


i(Wri)q(-w-) 


2/2 


CR2tS)  “  3tT3S2 


% 


d0  sin0  _ 
wl  wl 


r00 

„  2 

_  2 

,TT 

2  — Wp 

dw„w.  e  ^ 

dWlWl  6  ^ 

d8  _sin0  . 

2  2 

w2  w2 

1 

0 

TT 

f 

0 

2U 

0 

2TT 


(3-49c) 


d<fr  i  r| 

Wl  TT 


d0 *sin0 ’ IcosO * 
cr  ct‘  a 


lswuV^  <wR2)e 


% 


C3-49d) 


For  various  values  of  the  speed  ratio  S  the  collision  integrals  were 

+ 

again  evaluated  numerically.  The  results  obtained  by  averaging  the  inte¬ 
grand  over  40,000  random  points  are  presented  in  Tables  III  and  IV.  From 
the  figures  presented  in  Table  III  we  note  that  the  contributions  C  and 

xi2 

CR2  are  they  decrease  rapidly  with  increasing  values  of  the  speed 

ratio  S  and  can  be  neglected  for  S  >.  2.  The  results  for  S  =  0,1  are  in 
satisfactory  agreement  with  the  low  velocity  values  earlier  presented  in 
Table  II,  while  at  high  velocities  the  results  do  indeed  approach  the 
exact  values  given  in  C3-42) . 

^Further  details  concerning  the  procedures  used  in  calculating  the 
collision  integrals  cure  documented  in  a  Ph.D.  thesis  to  be  submitted  by 
Y.  Y.  Lin  Wang. 
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TABLE  III 

Collision  integrals  for  the  drag  coefficient  of  a  disc  (SSL2) 


s 

CH1XS 

CH2XS 

CR1XS 

C  XS  \ 

R2 

mm 

0.1 

-3.89±0.11 

+0.21±0v02 

+1.16±0.09 

-0.04+0.03 

-1.44±0.20 

0.5 

-4.96+0.15 

+0.16+0.02 

+1.42±0. 11 

-0.03+0.02 

-2.02±0.25 

1.0 

-8.01±0.16 

+0.068+0.005 

+2 . 16±0 . 10 

-0.014+0.007 

-3.65+0.22 

1.5 

-12.3+0.7 

+0.015+0.002 

+2.7±0 . 2 

-0 . 0019±0 . 0028 

-6 . 8±0 . 7 

2.0 

-15.0+1.4 

+0 . 0019+0 . 0005 

+3. 8+0. 6 

-0.0005+0.0003 

-7 . 4+1 . 6 

.  TABLE  IV 

Collision  integrals  for  the  drag  coefficient  C^  of  a  disc  (S>2) 
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3.7  Discussion  of  Results  for  the  Drag  of  a  Disc. 

In  summarizing  our  results  for  the  drag  in  the  nearly  free  molecular 
regime,  we  write  the  drag  coefficient  CD  as 


°d-  co[i+  5;  k'1]  ■  ^K^r8)*"1] ' 


(3-51) 


with  the  inverse  Knudsen  number  K  ^  defined  in  (3-10) .  The  term  (C,/C  )K  1 

1  o 

represents  the  first  inverse  Knudsen  number  correction  to  the  drag  force 
relative  to  the  magnitude  of  the  drag  force  in  the  free  molecular  flow 
limit.  The  contribution  C  =  C„.+C  accounts  for  the  effect  of  the 

ii  HI  HZ 

hypothetical  collisions  and  the  contribution  C  =  C  .+  C  for  the  effect 

R  R1  RZ 

of  the  recollisions  and  cyclic  collisions. 

We  have  calculated  the  coefficients  C  and  C,  over  the  entire  range 

o  1 

of  speed  ratios  between  zero  and  infinity,  assuming  that  the  temperature 
of  the  reflected  molecules  is  the  same  as  the  temperature  of  the  molecules 
in  the  incident  gas  stream.  The  free  molecular  flow  drag  coefficient  Cq 
is  given  by  (3-20) ,  while  the  coefficient  of  the  first  inverse  Knudsen 
correction  was  evaluated  numerically.  The  results  are  summarized  in 
Table  V  and  plotted  graphically  as  a  function  of  the  speed  ratio  S  in 
Figs.  10  and  11. 


In  the  nearly  free  molecular  flow  regime  the  drag  coefficient  decreases 

with  decreasing  Knudsen  number.  The  sign  of  the  effect  is  determined  by 

the  contribution  C  ,  i.e.  by  the  fact  that  the  reflected  molecules 

H  HI 

prevent  incident  molecules  from  striking  the  object. 

The  dependence  of  the  drag  coefficient  on  the  speed  ratio  is  quite 
different  whether  the  speed  ratio  is  smaller  or  larger  than  unity.  At  low 
velocities  the  drag  force  is  proportional  to  the  speed  ratio  S.  As  a 
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TABLE  V 


Drag  coefficient  of  a  disc  in  the  nearly  free  molecular  flow  regime  as 
a  function  of  the  speed  ratio  S. 


c 

o 

VCo 

C/C 

R  o 

c./c 

1  o 

0 

4.03 

-0.92±0.02 

+0.27±0.02 

-0.3810.02 

0.1 

4.04  S_T 

-0. 93+0.02 

+0.28±0.02 

-0.3610.05 

0.5 

4.21.S  J 

-1.1410.03 

+0.33±0.03 

-0.4810.06 

1.0 

4.72  S 

-1.6810.03 

+0.45±0.02 

-0.7710.05 

1.0 

4.72 

- (1.68±0.03)S 

+(0.4510.02)S 

-(0.77±0.05)S 

2.0 

3.14 

- (1.47+0. 01)S 

+(0.4410.01)S 

-(0.60±0.02)S 

3.5 

2.59 

-  (1.29±0.01)S 

+(0.4310.01)S 

-(0.43±0.01)S 

5.0 

2.39 

-(1.25±0.01)S 

+ (0.4310.01)S 

-(0.38±0.01)S 

7.5 

2.25 

- (1.19+0.01)S 

+  (0.4310.01)  S 

-(0.34±0.01)S 

10 

2.19 

-(1.16i0.01)S 

+(0.4210.01)S 

-(0.32±0.01)S 

15 

2.12 

- (1.14+0. 01) S 

+(0.4210.01)S 

-(0.29±0.01)S 

20 

2.09 

- (1.1110. 01) S 

+(0.4310.01)S 

-(0.26±0.01)S 

30 

2.06 

-(l.lOiO.Ol)S 

+ (0.4310.01)S 

-(0.26±0.01)S 

40 

2.04 

- (1.0910. 01)S 

+(0.4310.01)S 

-(0.23±0.01)S 

50 

2.04 

-(1.08i0.01)S 

+ (0.4210. 01)S  . 

-(0.23±0.01)S 

00 

2.00 

-1.064S 

+0r426S 

-0.213S 
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result  at  small  velocities  both  Cq  and  become  inverse  proportional 
to  the  speed  ratio  S  and  the  ratio  C^/Cq  becomes  independent  of  S.  On  the 
other  hand,  at  large  velocities  the  free  molecular  drag  force  varies  with 
the  square  of  the  stream  velocity,  while  the  first  inverse  Knudsen 
correction  to  the  drag  force  varies  with  the  third  power  of  the  stream 
velocity.  Thus  the  ratio  C^/CQ  becomes  proportional  to  the  speed  ratio 
S;  this  feature  is  encountered  for  a  large  class  of" objects  and  will  be 
further  discussed  in  Section  4.7.  In  the  limit  of  infinite  speed  ratio 
our  result  agrees  with  the  value  calculated  by  Willis  et  al.  [21]  for  the 
drag  coefficient  of  a  disc  in  a  beam  of  hard  spherical  molecules. 
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CHAPTER  IV 


DRAG  COEFFICIENT  OF  A  SPHERE  IN  THE  NEARLY  FREE  MOLECULAR  FLOW  REGIME 


4.1  Introduction 


As  a  second  application  of  our  method  we  consider  the  drag  force 
exerted  on  a  sphere.  A  number  of  investigators  have  attempted  to  calculate 
the  drag  coefficient  of  a  sphere  using  various  approximation  methods. 
However,  to  our  knowledge  a  solution  based  on  the  full  nonlinear 
Boltzmann  equation  that  covers  the  entire  range  of  speed  ratios  has  not 
been  presented  earlier. 

We  use  the  same  dimensionless  quantities  that  were  introduced  in 
Section  3.1  when  discussing  the  drag  coefficient  of  a  disc.  The  drag 
coefficient  in  the  nearly  free  molecular  flow  regime  is  again  represented 


by  (3-11) 


with 


C  =  C  + 
D  o 


C  X 
o 


-1 


-1  _  2 
K  =  \fl  TTnC  R 


(.4-1) 


(4-2) 


where  R  now  refers  to  the  radius  of  the  sphere.  The  drag  coefficient 
Cq  in  the  free  molecular  flow  limit  K  1-k>  is  in  analogy  co  (3-12) 


CQ  “  ITS' 


I  J*dw  f*  (w;S)  fdtk  |vft|  jdw'n*(wM  (w'-w) 


(4-3) 


w 


•ft<o 


The  coefficient  C^  of  the  first  inverse  Knudsen  number  correction 


V  °H1+  CH2+  2(CR14CR2)  ' 


(4-4) 
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is  again  determined  by  the  collision  integrals  listed  in  Table  I,  if 
?  and  n^  are  both  identified  with  the  unit  vector  in  the  direction 
of  and  and  with  the  unit  vector  &2  in  the  direction  of  R2  and 

where  R1  and  R2  are  indicated  in  Figs.  4  and  6.  The  dimensionless  velocities 

■+■  “V  *>■ 

Wi'  Wi'  Wi'  and  Wi''  ma^  with  the  velocities  v^,  v|,  v”  and 

vjl"  shown  in  Figs.  4  and  6,  provided  that  the  radius  of  the  sphere  is 
normalized  to  unity. 

4.2  Drag  Coefficient  CQ  of  a  Sphere  in  the  Free  Molecular  Flow  Regime. 

The  drag  coefficient  C  of  a  sphere  in  the  free  molecular  flow 

o 

A 

regime  is  given  by  (4*3)  •  Lets  be  the  polar  angle  of  R  with  respect 

R 

to  -S;  because  of  symmetry  we  can  integrate  over  the  azimuthal  angle 

A  A 

of  R.  Let  0^  be  the  polar  angle  of  w  with  respect  to  R  and  the 
azimuthal  angle  of  w  with  the  plane  through  S  and  R  as  initial  plane. 

Then,  using  (3-14) , 


C  (S). 
O 


/r  00  TT  27T 

d0Rsin0R/dw  w3/ d8wSin0W  lcosBJ  fdK 


lz  6 

l/rr 

(4-5) 

A  K 

with 

■+■  A 

w*S=  w(sin©  sin0 'cos^'-cos®  cos0')  .  (4-5a) 

R  w  w  R  w 

The  integral  can  be  evaluated  analytically  with  the  result  [1,10,32] 


.  e-(w2«2)+  2S<;-S>  {(J.g)+ScoaeR}  _ 


(4-6) 


In  the  low  and  high  speed  limits  the  drag  coefficient  C  (S)  reduces  to 

o 


lim  C  (S)  = 

s-+o  ° 


2TT+16  ^  4.191 
3\/f  S  S 


(4-7) 
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and 

lim  C  (S)=  2.  C4-8) 

S-x»  o 


The  coefficient  of  the  first  inverse  Knudsen  number  correction 

to  the  drag  coefficient  is  again  determined  by  the  collision  integrals 

r^- 

in  Table  I  with  r^=  n^=  and  n^-  R2*  T^e  that  each  collision 

integral  is  related  to  a  particular  collision  sequence  is  accounted  for 
by  the  d-functions.  In  this  Section  we  integrate  over  the  5- functions 
and  formulate  explicitly  the  conditions  for  the  occurrance  of  these 
collision  sequences. 

As  an  example  we  consider  the  collision  integrals  associated  with  the 
hypothetical  collisions  (cf .  Fig.  6) .  The  geometry  of  an  Hl-sequence 
is  indicated  in  Fig.  12.  The  Z-axis  is  taken  in  the  direction  of  -S  and 

A 

the  X-axis  in  the  plane  through  S  and  R^  Molecule  1  is  emitted  from  the 
surface  at  with  velocity  w^.  Molecule  2  with  incident  velocity  w2 
collides  with  molecule  1  after  a  time  T*  .  We  introduce  the  vector  r 
defined  as 

r  =  (4-9) 

which  determines  the  location  of  the  collision  between  molecules  1  and  2 


relative  to  the  center  of  the  sphere.  Molecule  2  continues  to  proceed 

with  its  initial  velocity  w  and  strikes  the  surface  of  the  object  after 

a  time  T*  .  Although  not  indicated  explicitly  in  Fig.  12,  the  vectors 

r  and  R2  are  not  restricted  to  the  XZ-plane. 

The  6-function  in  the  expression  (1-1)  for  C  can  be  rewritten  in 

HI 

-4  4 

terms  of  the  velocity  w2  and  the  distance  vector  r 


^Sections  4.3  and  4.4  were  prepared  in  collaboration  with  Dr.  W.  A.  Kuperman, 
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(4-10) 


(T2+ 


2  " 

r  -R 


¥■) 


6z(k-r-&Jw- 
2  2  2 


r+  t/l-  (w^r)  2})0(-w2*r)0(l-{w2*r}2)  . 


The  conditions  imposed  by  the  Heaviside  functions  in  (4-10)  have  an 

* 

obvious  geometrical  meaning.  The  condition  w2*x  <  0  guarantees  that 

•+■  2 

molecule  2  approaches  the  sphere  and  the  condition  (w2  x  r)  <1  ensures 
that  the  impact  parameter  is  smaller  than  unity, 

The  conditions  for  the  H2-sequence  (cf.  Fig.  6b)  can  be  related  to 
the  conditions  for  the  Hl-sequence,  if  we  identify  the  velocity  in 

*f 

(4-9)  and  (4-10)  with  the  initial  velocity  w^  and  if  we  register  the  first 
collision  not  at  the  position  where  molecule  1  enters  the  sphere,  but  at 
the  position  where  molecule  1  leaves  the  sphere. 

^  A 

Hence,  if  we  integrate  over  T*  ,  R2  and  CT12(cf. (3-32) )  and  use  (3-14), 
the  collision  integrals  (1-1)  and  (1-2)  reduce  to 


CH1(S)=  -jp|j  Jdw1dl2f*(w1;S)f»(»2iS)  JaT*  JdHj  |w1'91l 


A 

W1 *R1<0 

dw|Tl*  (wpw|2  (w2-  2).g  0(-w2-?)0(l-  {w2*r}2) 


(4-lla) 


CH2  (s)  =  +  H2  / dV«2f*  f*  ( fc*l  / ®1  I 


w  .ft  >° 


’W12<”2_J^‘  ^2J  0("W2*r)0(l-{w2><r}2) , 


(4-llb) 


with 


ft. 


r  r_w2|w; 


•r+ 


»2Xr)' 


) 


(4-llc) 


The  collision  integrals  (1-3)  and  (1-4)  can  be  treated  in  the  same 

manner,  except  that  we  cannot  integrate  analytically  over  the  collision 

vector  d. Thus 
12 
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oo 

®R1CS)=  +tSp  tWi;S)f*  (w2;S)  f  dT*  /dR1  |w 

w^ • R^<o 

•Jdwjn*(w^)  j  d&12  |w^2-ei2|  Cw|‘-^~ft2).§  0  (-^' *r)0(l-{w^ 'Hr}2) ,  (4- 

"i2*12<0 

CR2(S)=  "  ^S2  fryS)  f&2;S)  fitt*  fd^ 

•  n  *  4-  »  . 


w^«R^>o 


■/d912  I”l2‘d12l  ^  &2)  -S  Sl-W^'  .  (4- 12b) 


W12 *^12<0 


with  .  ^ 

&2=  r-w^,\w|,,r4  l/l-  (w|'Xrj2  /  . 


(4-12c) 


The  velocity  w|‘  is  again  related  to  the  integration  variables  by  C3-35) . 


4.4  Drag  Coefficient  of  a  Sphere  at  High  Speed  Ratios. 

In  the  high  speed  limit  or  (cold  wall)  beam  limit 

f* (wi;S)=  6 3 (w^-S)  ,  (4-13) 

and 

lim  C. (S)=  lim  c  (s)+  2  lim  C_n  (S)  ,  (4-14) 

S^oo  1  g-*»  HI  s-*50  R1 

as  discussed  in  Section  3.4.  Retaining  only  the  leading  term,  in  S  the 
collision  'integrals  (4-lla)  and  (4-12a)  reduce  to 


lim  C 

Sx» 


HI 


f dT*  fd^  IS-fcJ  ydw|n*(wp0(-§-r)0(l-{§Xr}2),  (4-15) 
0  S  •  Rx<  o 


lim  C 
S+“ 


R1 


OO 

(S)= +§s  Jdt*  18^1  J^n*Cw|) J d&l2(S  di2) 

§*ft  <o  S*d  >o 

1  12 

•G(-ei2^)0(l-{&12x^}2)  . 


3 


(4-16) 


12a) 
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We  express  the  vectors  ^(6^4^) ,  r(r r0r ,4*^)  and  »4>wli J 

in  terms  of  spherical  polar  coordinates  in  the  coordinate  system  with 
the  Z-axis  in  the  direction  of  -S.  Because  of  symmetry  we  may  integrate 

.  A 

over  the  azimuthal  angle  <J>R1  of  and  take  in  the  XZ-plane  as  indicated 
in  Fig.  12.  Then 

a>  jsfo  00  ,2  IT  2Tr 

Hl(S)=“f  S  JdTlJd0RlSin9RlC°S0Rl  JdWiWi3e’Wl  \  d0wl'Sin0wH  jdKv 

0  0  0  0  0 

*(w' ’ft  )0(w'  )0(cos9  )0 (1-rsinS  )  ,  (4-17) 

x  x  i  l  r  r 

with 

cos0R1cos0wlI+sir6R1sin9wl,cos(()wll  .  (4-17a) 

The  Heaviside  functions  in  (4-17)  impose  the  conditions 


°<0R1<  arct?T  t-cotg0wl,  sec<t>wi,  )  »  (4-18) 

and 

0.<©r  <arcsln  r  1 .  C4-19) 

It  follows  from  (4-9)  that  condition  (4-19)  is  equivalent  with 


0<T1*<  81110  Hi 


w'sin0  , 
1  wl 


■  [ 


-  c°a4>  wl-+  ^  co  tg  26  Ri+  cos2<Pwl- 


I  • 

•  • 


(4-20) 


With  the  integration  limits  for  T*  and  0  determined  by  (4-20)  and 
(4-18)  the  Heaviside  functions  in  (4-17)  are  taken  to  be  unity.  All 
integrations  can  be  performed  analytically  with  the  result 

«-21> 


In  order  to  evaluate  the  collision  integral  (4-16) ,  as  well  as  the 
collision  integrals  in  the  subsequent  sections,  we  consider  three 
different  coordinate  systems,  to  be  referred  to  as  the  XYZ,  X'Y'Z’  and 
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X* • Y*  *Z* '  systems.  The  XYZ  system  is  the  original  coordinate  system 
with  the  Z-axis  in  the  direction  of  The  X'Y'Z'  system  has  the  Z'-axis 

in  the  direction  of  1?  and  the  X'-axis  in  the  plane  through  and  §.  The 
X' ' Y'  'Z * '  system  has  the  Z"-axis  in  the  direction  of  r  and  the  X"-sxis 
in  the  plane  through  r  and  S.  We  express  w{  (wi'®wl' '^wl1 >  -in  terms  of 
polar  coordinates  in  the  X’Y'Z1  system.  The  polar  and  azimuthal  angle  of 
in  t^ie  original  system  are  0ff,  and  in  the  x^Y^Z'*  system 

■f 

The  Cartesian  components  of  any  vector  a  in  the  auxiliary  systems  are 
related  to  its  components  in  the  original  XYZ  system  by 


/a  \ 

fa  \ 

/a  ^ 

fa  \ 

X' 

x\ 

ax"’ 

ax\ 

V 

Ry(0Ri} 

aY  ' 

ay" 

-  W  W 

a 

y 

\'"l 

W 

f-7 

laJ 

(4-22) 


where  R  and  R  are  the  same  rotation  matrices  as  those  used  in  earlier 
j  z 

reports  [3i,33].-and  where  0^  and  0^  are  the  polar  and  azimuthal  angle  of 


the  distance  vector  r 


r  =  R  +  w't*=  ft  +w 
1  11  1  1 


'  f-cos0V 
1  I  w. 


,+s/r5- 


sin^0 ' 


,  i  ■  v  j.  v  -  I 

wl*  wl' 

The  collision  integral  (4-16)  reads  in  terms  of  these  variables 

00  m 


(4-23) 


lim  C 
3*“  R1 


(S)=  +  *^sjax*  Jde^sinO^cose^  jdw- w- 3e_wi  2  J?0^ifsiB0;i,cos0;1( 
2TT  m.  ZTT 

•^d(|>,wlIJd0CTs^h0acos30o.^|,d<^)(70(-Qos0(:J, )  0(l-rsinl 


rsin0 ' ' ) . 

o 


The  Heaviside  functions  in  (4-24)  impose  the  condition 

ir-  arcsin  r  1  <  0”  <  n  ,  (4-25) 

so  that  molecule  1  will  indeed  return  to  the  sphere  after  the  collision 
with  molecule  2.  We  transform  the  integral  over  the  time  T*  from  0  to 
00  into  an  integral  over  the  distance  r  from  1  to®  [19] 
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dr*= 


rdr 


1  Vr^-sin2©^, 


(4-26) 


If  we  then  integrate  over  w^,  we  obtain 


lint  C  (S)=  + 

S-x»  R1 


'-i.s  f  ae„,sin9„,oose 


■jdra-r-2sin2e^-1/2 


Rl  UR1 

4 


d9  sin0  cos  0_ 

a  a  a 


Rl/ kv 


2TT 


sin0 ' , . cos0  , 
wl'  wl 


ir  /  **h- 


2-T 


0(-cos0',)&(l-rsin0")  .  (4-27) 
O  0  o 


The  remaining  integral  was  estimated  numerically  as  a  weighted  average 
over  100,000  random  points  with  the  result  +(0. 510*0. 002)S. 

We  thus  conclude  that  in  the  limit  of  high  speed  ratios  [16] 


&  CB1<S)“  'Vi  S  =  -1'253  S  < 

Jim  CRl(S)=  +(0.510±0.002)S  , 

J^ffiC1(S)=  Jim  fcH1(S)+2CRl(S)}  -  -(0.233±0.004)S  . 


(4- 28a) 
(4- 28b) 
(4- 28c) 


These  results  are  in  good  agreement  with  the  values  Jdgj  CH1(S)  = 

- (1.23+0. 02)S  and  Jim  CRl(S)=  +(0.511+0.002) S  as  calculated  by  Kupermant. 


t  Equation  (.5-54)  in  the  original  thesis  of  Kuperman  contains  an  error  and 

should  read  v^aj™  ^x^^z^wl^Oci)  *  Th*s  chan9e  does  not  affect  the 
value  of  ClHl  quoted  in  Table  5.2  of  his  thesis,  but  the  values  quoted  for 
CXR1  and  C^1  should  be  replaced  with  C1B1^1C1=(.0,256±0.002IS/Kn' [35J  ,  For 
the  same  reason  the  values  quoted  in  Table  6.4  of  Kuperman’s  thesis  should 
be  replaced  with  the  values  presented  in  Section  4,5  of  this  report. 
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4.5  Drag  Coefficient  of  a  Sphere  at  Low  Speed  Ratios. 

At  low  speed  ratios  we  can  expand  the  Maxwell  distributions  in  the 
collision  integrals  (4-11)  and  (4-12)  in  a  Taylor  series  around  S=0  and 

retain  only  the  first  order  term  in  S. 

+  +  +  +  1  —  (w?+w2)  r  , 

f*(w  »S)f*(w  ;S)=  -?  e  '  1  2;{1+2S.(W  +w  )}  .  (4-29) 

1  2  7r3  12 

The  zeroth  order  term  in  (4-29)  does  not  contribute  to  the  drag  coefficient. 

We  use  the  coordinate  systems  introduced  in  the  preceding  section. 

The  polar  and  azimuthal  angle  of  any  vector  a  are  denoted  by  0  ,<J>  ,  0 1  ,tf» 1 

d  d  cl  cl 

and  in  the  XYZ,  X'Y'Z*  and  X^Y^Z11  system,  respectively.  It 

a  a 

is  also  convenient  to  introduce  the  polar  and  azimuthal  angles  0^  and  <|>a  of 
in  a  coordinate  system  with  the  Z-axis  in  the  direction  of  w'  in 
(4-12a)  and  in  the  direction  of  w12  in  (4-12b)  and  the  X^-axis  in  the  plane 
through  S  and  w|2  or  w12-  We  transform  again  the  variable  T*  into  the 
distance  r  in  accordance  with  (4-26) .  In  the  collision  integral  (4-lla)  we 
can  integrate  readily  over  the  velocity  w^  We  thus  obtain 

2  00  2  IT 

/of  2  — wj *  r 

C„,  (S)=  -  -r:  | dw„w~e  2  I  dw>^  e 

&*-  o 


i 

/ 


CHllS)*  -  ^/d”2U2e"“2/d"iWi2e‘Wi  /d3RlSl"eR1 

0  0  0 

J/2  2lT  00  TT  2ir 

d0wl'sin6ul'cosewl'  f^wl'  /ara-r‘2sinVlirl/2  Ja9-Sin8-  Jd*- 

o  i  - J-“-1  0 

‘f  R2z'”2z,(  J?“"6KL'*te,r 

c 

■r< 


tf-arcsinr 


W12< 


(4-30a) 


4  /2 

lun  C  (S)=  4^ 
S+  9  H2  TT sj 


00  _  00 

.  2  -wf 

dWlWlS 


JdW2V2e~W2  f 

0  0 


d0  sin8 
R1  Rl 


•IdS'  sin0'  cos0  , 
J  wl  wl  wl 


2ir 


il/d*wl/dr<1"r 


-r_2sin20 ') "1/2 


7T 


2TT 


ii rl/2  f* 

TT-arcsinr”  r 


I  I 

w2 


•»12<''22-^|'>2z)(«i2+»2z). 


(4-30b) 
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lilt*  C„.  (S)  = 


oo  2  co  2  7T  7T  2.7T 

’’  +  4?/aV2e'“2/aWi"i2e'Wi  /d0RISln6Rl/d02Sirie2/<i+2 
0  o  0  0  0 

TT^  2TT  a,  IT  *7T 

'j^e*l,sinewlcos0wl,  /^il,  f  ar(1'r’2«in29^ll>'1'!,2/a5osiii0a  |=os8  |  (ai 
0  0  1  %  o' 

'Wl2(-VWu)  (_T  cos9r1-V°s62)  e("COS0il,r)0<l-rsin0^II)  ,  (4-30c) 


lim  C__ 

S*0  R2 


_  00  2  •  2  IT  7f  21T 

<S>-  -  ^/nV^1j*'2”2e'W2/dSR1Sl"eRldfde2si,«2/'il'2 


M  27T  00  7T  air 

J®  wi3irfi  liooEe  ;,i  /dr  (1-r_2Bi'A  ;*>  "1/2  / *  |™=80I  Ua 

0  0  1  7%  O 


•wl2(wiz‘  i?R2z)  (wlZ+W2z)®  {-COsewi  *  •^(l-rsin0^, , )  . 


C4-30d) 


These  integrals  were  again  calculated  numerically.  The  results  are 
presented,  in  Table  VI. 


4,6  Drag  Coefficient  of  a  Sphere  at  Arbitrary  Speed  Ratios 

For  arbitrary  values  of  the  speed  ratio  S  we  need  to  substitute  the 
complete  Maxwell  distribution  (3-6)  into  the  collision  integrals  (4-11)  and 
(4-12) .  In  formulating  the  collision  integrals  we  use  the  same  notation  as 
in  the  preceding  section  and  obtain 


„  4/2 

CH1(S) 


,TT 


r°° 

.00 

f<x> 

2 

„  2 

n  ,2 

.  3  -wi 

dw,  w,  e  A 

2  -w, 
dw_w_  e  2 

.  .  ,2  -wl 
dw'wl  e  1 

1  1 

2  2 

1  1 

oJ 

0^ 

oJ 

r  \v> 

,2Tr  ,00 

f  r 

r* 

d0RlSin0Rl 


d0'sin0'  cos0' 
wl  wl  wl 


% 


dKl 


d9Wl'Sill0«7  ■COS0-T 


wl1 


wl' 


'wl 1 


dr ( 1-r” 2  s in20 '  ,)“1/2 
wl 


rl T 

d0 ' ' sin0 ' ' 
w2  w2 


^..e-2s<s«,1^,w.2  , 


TT-arcsin  r 


-1 


(4-31a) 
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TABLE  VI 


Collision  integrals  for  the  drag  coefficient  of  a  sphere  in  the  low 
velocity  limit. 


lim  C. (S)=  -(2.585±0.014)S 
„ .  HI 

s+o 

lim  C  _(S)  =  +(0.299±0.008)S-1 
S+O 

lim  CH(S)=  (3)^(3)} 

CR1(S)=  +  (0.59±0.03)S~1 

CR2(S)=  - (0.059-0.007) S_1 

lim  c  (s)  =  lim  {c_.  (s)+c_.  (S) } 
s-K)  R  s+o  R1  R2 


=  -(2.29±0.01)S 


=  +C0.53±0.03)S 


lim  C. (S>=  lim{C  „(S)  +  2C  (S)}=  - (1.23+0.05) S 

s+o  1  s+o  H  R 
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CH2(S)  *  +  510 


-  2  -w, 

dwiwi  e 


,  2  -w, 

dw2w2  e  * 


7T 

deBlSin9Rl 


yi 

2ir 

#00 

3 

• 

d0 1  sin0 ' .  cos0 • . 
wl  wl  wl 

d(j> ' 
wl 

,  -2  .  2ot  .  -1/2 
dr (l-r  sin  0'  ) 

Wi 

dO'lsinO': 
w2  w2 

0 

A 

iJ 

-1  nJ 

2ir 


d^'l 
Yw2 


ir-arcsm  r 


-2S (S+w 


’  +w  )  //F  _  \ 

12  22  W12  It  R2z  -  W2zJ  ■ 


(4-3IM 


4^2" 

CR1(S>  “  +  TTsS2 


O 


.00 

fOO 

_  2 

_  2 

„  ,2 

dw  w  e 

dw  w  e  W2 

,  ,2  -w' 
dWll  6  3 

1  1 

2  2 

TT 


d0;isin0;iicos0;ii 


% 


27T  .TT 


d02sin02 


o  0 

rir 


0 

air  % 


„  „  -2  .  2Qi  .-1/2 
dr  (l-r  sin  0^,) 


1' 


d*2 

o'  o' 


TT 


d0iasineKL 


2ir 


d0 '  , sin0 ' n , cos0 ' .  , 
wl 1  wl '  wl ' 


d0CTsin0a| cos0a| 

74  0 

•  0(-cos0^I1)0{l-rsin0^|ll) , 


~  “2S (S+w.  +w  ,,\ 

“V  iz  2z  R22-»i;J 


(4-31c) 


C  (s)  _ ?/[2 

CR2l  J  TT**S^ 


.  2  -w-i 

dwiwi  e  1 


dw2w22e_W2 


d0Rlsin0Rl 


oJ 

oJ 

1 

2-rr 

.IT 

2TT 

.CO 

•  d0 '  sin0 '  cos0 ' 
wl  wl  wl 

d(j>' 

wl 

d02sin02 

d<j>  ~ 
rw2 

.  -2  .  20  -1/2 
dr (l-r  sin  0*  ) 
wl 

TT 


d0  sin0„,  cos0„ 

a  o'  o 


0  0 

'27r 

d<Le_2S(S+Wlz+W2z)w,  R„  -w„"W- 


* 


12\  2  *'2z  "izj0(_cos9wi")0(1-rsJ'n0wr,) 


(4-3 Id) 


The  various  integrals  were  again  computed  numerically.  In  Table  VII  we 
present  the  estimated  values  of  C^,  as  well  as  the  individual  contributions 

CHlf  CH2'  CRlf  CR2  for  S-3,  Just  as  in  Section  3.6  when  calculating  the  drag 
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TABLE  VII 


Collision  integrals  for  the  drag  coefficient  C^  of  a  sphere  (S<3) 


s 

C  .*S 
«1 

c  „xs 

H2 

CR1XS 

C  *S 

R2 

c^s 

0.01 

-1.22+0.10 

0.05 

-2.51+0.07 

+0. 26±0.07 

+0.4510.03 

-0.0910.03 

-1.1610.06 

0.10 

-1.2410.10 

0.50 

-2 . 94±0 . 16 

+0.3010.09 

+0.6510.05 

-0.07+0.02 

-1.4710.16 

1.00 

-4. 51±0.12 

+0.1410.01 

+1.0410.05 

-0.045+0.006 

-2.3610.14 

1.50 

-6.47±0.25 

+0.05410.005 

+1.48+0.07 

-0.016+0.002 

-3.4810.27 

3.00 

-17.95+0.39 

+0.002+0.001 

+5.4710.39 

-0.0008+0.0004 

-7.0010.85 

TABLE  VIII 


Collision  integrals  for  the  drag  coefficient  C.^  of  a  sphere  (S£5) . 


S 

CH1/S 

CR1/S 

c1/s 

5.0 

-1.6410.04 

+0.5210.02 

-0.6010.05 

7.5 

-1.4910.03 

+0.5110.01 

-0.4910.04 

10. 

-1.44+0.03 

+0.5010.01 

-0.43+0.03 

15. 

-1.3210.03 

+0.52+0.02 

-0.3710.04 

20. 

-1.2910.03 

+0.5210.02 

-0.2510.04 

30. 

-1.26+0.02 

+0.49+0.01 

-0.29+0.03 

50. 

-1.26+0.05 

+0.50+0.01 

-0.2610.03 
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of  a  disc,  we  again  notice  that  the  contributions  c_  and  C__  may  be  neglected 
for  speed  ratios  S>2.  The  results  obtained  for  the  larger  speed  ratios 
are  presented  in  Table  VIII.  At  low  and  high  speed  ratios  the  results 
approach  the  limiting  values  established  earlier  in  Sections  4.4  and  4.5. 


4.7  Discussion  of  Results  for  the  Drag  of  a  Sphere. 

In  summarizing  our  results  for  the  drag  of  a  sphere  in  the  nearly  free 
molecular  flow  regime,  we  write  the  drag  coefficient  CQ  again  as  in  (3-51) 


C  =  C 
D  o 


[■•ri-M-ffr] 


(4-32) 


with 


K_1=  /2  irncr2R. 


(4-33) 


A  summary  of  the  values  of  the  coefficients  in  the  expansion  (4-32)  for  the 
drag  coefficient  CD  is  presented  in  Table  IX.  The  coefficients  Cq  and 
C^/C  are  plotted  as  a  function  of  the  speed  ratio  S  in  Figs.  13  and  14 
using  the  same  scale  as  used  earlier  for  the  drag  coefficient  of  a  disc  in 
Figs.  10  and  11.  On  comparing  Table  IX  with  Table  V  and  Eigs.  13  and  14 
with  Figs.  10  and  11,  we  note  that  the  drag  coefficient  of  a  disc  and  a 
sphere  are  rather  similar  functions  of  the  speed  ratio.  Again  at  low  speed 
ratios  C^/Cq  becomes  independent  of  the  speed  ratio,  while  at  large  speed 

ratios  C,/C  becomes  proportional  to  the  speed  ratio  S. 

1  o 

The  drag  exerted  on  a  sphere  in  the  nearly  free  molecular  flow  regime 
has  been  studied  by  a  number  of  authors.  Host  of  these  studies  are  concerned 
either  with  low  velocities  S<<1  or  large  velocities  S»1  and  we  shall  discuss 
the  two  cases  separately. 

Our  results  may  be  interpreted  as  the  solution  of  the  Boltzmann  equation 
for  a  gas  of  hard  spheres  in  the  presence  of  the  object.  A  study  of  the 
sphere  drag  based  on  the  Boltzmann  equation  for  Maxwellian  molecules,  i.e. 
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TABLE  IX 


Drag  coefficient  of  a  sphere  in  the  nearly  free  molecular  flow 


regime  as  a  function  of  the  speed  ratio  S. 


C 

o 


k_1=  V2Trna2R. 


s 

Co 

Cr/Co 

<Vco 

ci/c0 

0 

4 . 19S^f 

-0. 546±0. 003 

+0.12610.006 

-0.2910.01 

0.01 

4.19s'1 

-0.527±0.007 

+0.1210.01 

-0.2910.02 

0.05 

4.19S_; 

-0.539+0.003 

+0.1310.01 

-0.2810.01 

0.10 

4.20S_1 

-0. 52910.008 

+0.1310.01 

-0.2810.01 

0.50 

4.34S  1 

-0.61±0.03 

+0.1310.01 

-0.3410.04 

1.0 

4.75S  1 

-0 . 9210 .02 

+0.2110.01 

-0.5010.03 

1.0 

4.75 

-  (0.92±0.02)S 

+ (0. 2110.01) S 

-(0.5010.03)S 

1.5 

3.58 

-(0.80±0.03)S 

+ (0.18±0 .01) S 

-(0.4310.03) S 

3.0 

2.61 

-  (0.76+0. 02)  S 

+  (0.23+0.02) S 

- (0.3010.04) S 

5.0 

2.32 

- (0.71±0.02)S 

+T0.2210.0DS 

- (0 . 2610.03) S 

7.5 

2.19 

- (0.68+0.01) S 

+(0. 23010. 005)S 

- (0.2210.02) S 

10. 

2.14 

-(0.67±0.01)S 

+(0. 23510. 006)S 

-(0.2010.02) S 

15. 

2.09 

-(0.63±0.01)S 

+(0. 24810.008) S 

- (0.1310.02) S 

20. 

2.06 

- (0.63±0.01)  S 

+(0. 253+0. 008)S 

-(0. 12+0.02) S 

30. 

2.04 

- (0. 63±0 . 01)  S 

+(0.24+0.01)5 

-(0.1410.01)  S 

50. 

2.02 

-(0.6210. 02)S 

+ (0.24710. 008) S 

-(0.1310.01)  S 

00 

2.00 

-0.627S 

+(0. 25510. 001)S 

-  (0.11710. 002)  S 
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Figure  14,  The  coefficient  Cj/Cq  of  the  first  inverse  Knudsen  number 

correction  to  the  drag  coefficient  of  a  sphere  as  a  function 
of  the  speed  ratio  S. 
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molecules  that  repel  each  other  with  forces  proportional  to  the  inverse 
fifth  power  of  the  intermolecular  separation,  was  made  by  Liu,  Pang  and 
Jew  [4] .  They  expanded  the  collision  integrals  into  Hermite  polynomials 
in  terms  of  the  molecular  velocity;  this  procedure  leads  to  an  expansion 
of  the  drag  force  in  which  higher  order  terms  of  the  speed  ratio  S  are 
neglected . 

A  comparison  of  our  results  with  those  of  Liu  et  al.  is  made  in  Table  X. 
The  predicted  drag  coefficients  appear  to  be  very  similar  for  S<0.5;  in  this 
range  we  confirm  the  conclusion  of  Liu  et  al.  that  the  coefficient  C]/CD  is 
almost  independent  of  the  speed  ratio  S.  For  £>>0.5  the  results  of  Liu  et  al. 
begin  to  deviate  substantially  from  our  results;  this  difference  might  be 
due  to  a  failure  of  the  expansion  procedure  of  Liu  et  al.  for  larger  values 
of  £>  [4].  The  convergence  of  the  expansion  procedure  of  Liu  et  al.  has  been 
questioned  by  Willis  [7].  Nevertheless,  unless  the  agreement  is  fortuitous, 
the  data  of  Table  X  would  suggest  that  the  drag  coefficient  is  insensitive 
to  the  details  of  the  molecular  interaction  and  that  it  is  mainly  a  function 
of  the  size  of  the  interaction  range,  or,  equivalently,  of  the  magnitude 
of  the  mean  free  path. 

Due  to  the  complications  associated  with  solving  the  full  Boltzmann 
equation  many  authors  have  used  instead  the  BGK  equation  (1-4)  or  more 
sophisticated  versions  of  this  model  equation.  The  use  of  this  model 
equation  Introduces  an  uncontrolled  approximation  for  mathematical  conve¬ 
nience.  Its  predictive  power  is  further  limited  by  the  appearance  of  an 
adjustable  collision  frequency  V. 

In  the  limit  9*o  it  is  sufficient  to  consider  a  linearized  version  of 
the  BGK  equation  (1-4) .  Using  this  procedure  and  assuming  again  diffusive 
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TABLE  X 


Comparison  of  our  results  for  the  sphere  drag  with  those  of  Liu  et  al. 


Liu  et  al  [4] 
(Maxwellian  molecules) 

This  work 
(hard  spheres) 

S=  0 

C-/C  =  -0.298 

1  o 

C,/C  =  -0.29±0.01 

1  o 

0.01 

-0.298 

-0.29±0.02 

0.10 

-0.300 

-0.28±0.01 

0.50 

-0.308 

-0.34-0.04 

1.00 

-0..296 

-0.50±0.03 

TABLE  XI 


Survey  of  theoretical  values  reported  for  lim  C../C  of  a  sphere. 

1  ° 


Cl 

lim  —  =  -0.24  S 
S-w  C 

o 

=  -0.143  S 


=  -0.33 


=  -0.165 


r  Baker  and  Charwat  [14,38] 
,  Perepukhov  [15] 

/  Rose  [8,39] 

/  Willis  [38] 


=  -(0.117± 0.002)  S  ,  this  work. 
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reflection  of  the  molecules  from  the  object,  Willis  calculated  the  coefficient 
C^/Cq  for  a  sphere  with  the  result  [7J 

lira  Si  =  -0.366.  (4-34a) 

S*0  Co 

Starting  from  the  same  equation  we  find^ 

lim  =  - (0.4Q1±0.004)  .  (4-34b) 

S+O  C 

o 

These  theoretical  values  are  usually  compared  with  the  measurements  of 
Millikan  [36]  for  the  drag  of  oil  droplets  in  air  as  a  function  of  the 
Knudsen  number.  From  the  experimental  data  of  Millikan  one  may  deduce^ 

C1 

0.39+0.02  -  (4-35) 

Co 

The  rather  close  agreement  between  (4-34)  and  (4-35)  has  been  widely  in¬ 
terpreted  as  a  justification  for  the  use  of  the  BGK  equation.  However,  such 
a  conclusion  is  dangerous,  since  it  heavily  depends  on  the  presupposition 
that  the  theoretical  boundary  conditions  are  satisfied  in  Millikan's  ex¬ 
periment.  In  fact,  Millikan' s  experiments  were  conducted  in  air  which  is  a 
mixture,  while  the  results  are  interpreted  in  terms  of  equations  for  a  one 
component  system.  A  more  reliable  criterion  for  the  adequacy  of  the  BGK 
equation  is  obtained  by  comparing  its  solution  with  that  of  the  Boltzmann 


This  work  will  be  included  in  a  Ph.D.  thesis  to  be  submitted  by 
Y.  Y.  Lin  Wang. 

■ft 

A  statistical  analysis  of  Millikan's  data  in  terms  of  the  theoretically 
predicted  equation  CQ=  was  ma^e  ^  Zondlo  [373  . 
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equation  with  the  same  boundary  conditions  [12] .  The  fact  that  the  BGK 
result  (4-35)  differs  appreciably  from  the  Boltzmann  equation  result  both 
for  hard  spheres  and  Maxwellian  molecules  raises  doubt  about  the  adequacy 
of  the  BGK  model  equation. 

Theoretical  calculations  of  the  sphere  drag  at  large  velocities  have 
been  presented  by  Baker  and  Charwat  [14],  by  Perepukhov  [15],  by  Willis  and 
coworkers  [38]  and  by  Rose  [8,39] .  A  summary  of  the  values  reported  for 
^im  Ci/CQ  is  presented  in  Table  XI.  In  order  to  simplify  a  comparison  with 
our  results  we  have  only  quoted  in  Table  XI  the  limiting  values,  when  the 
temperature  of  the  object  is  identified  with  the  temperature  of  the  gas 
stream.  Baker  and  Charwat  [14]  and  Perepukhov  [15]  start  from  the  same  model 
as  considered  in  this  report,  but  then  proceed  to  introduce  a  number  of 
approximations.  The  extent  to  which  their  results  differ  from  our  value  of 
- (0.117 ±0.002) S  indicates  the  effect  of  their  approximations.  In  particular 
the  drastic  approximations  introduced  by  Baker  and  Charwat,  do  not  appear 
to  be  justified.  The  fact  that  Cj/C0  becomes  proportional  to  the  speed  ratio 
S  at  high  velocity,  is  a  consequence  of  the  general  nature  of  the  collision 
integrals . 

The  value  obtained  by  Rose  [8,39]  is  based  on  the  BGK  equation  (1-4) 
and  that  obtained  by  Willis  [38]  is  based  on  a  modified  version  of  the  BGK 
equation.  The  dependence  of  their  solutions  on  the  speed  ratio  S  appears 
to  be  qualitatively  different.  These  model  equation  results  are  subject  to 
the  limitations  mentioned  earlier  in  the  discussion  of  the  drag  coefficient 
in  the  low  velocity  limit.  Their  validity  can  only  be  judged  a  posteriori 
from  a  comparison  with  the  solution  of  the  Boltzmann  equation,  though  not 
necessarily  the  solution  of  the  Boltzmann  equation  for  hard  sphere  molecules. 
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The  fact  that  our  calculated  values  of  C]/CQ  become  proportional  to  the 

speed  ratio  S  at  large  velocities  deserves  some  further  comments.  The 

validity  of  the  expansion  (4-32)  for  the  drag  coefficient  is  limited  to  the 

range  where  the  mean  free  path  is  substantially  larger  than  the  object,  i.e. 

the  inverse  Knudsen  number  should  be  substantially  smaller  than  unity.  How- 

—12 

ever,  the  parameter  K  =  V2  'UnO  R  in  (4-32)  represents  the  Knudsen  number  in 
the  absence  of  (i.e.  away  from)  the  sphere  and  is  therefore  sometimes  referred 
to  as  Km1,  in  order  for  the  expansion  to  be  valid  we  must  require  that  the 
local  inverse  Knudsen  number  near  the  object  is  substantially  smaller  than  .. 

•v 

unity.  Indicating  the  local  Knudsen  number  by  K,  then  at  small  velocities 

lim  "k_1=  K_1=^  irrxr2R  ,  (4-36a) 

9+o 

while  at  large  velocities 

lim  K-1^  SK_1=  Si/T 7Tna2R  ,  (4-36b) 

S-Xn 

since,  at  large  velocities,  the  mean  free  path  of  the  reflected  molecules 
becomes  inversely  proportional  to  the  speed  ratio  S.  If  we  thus  rewrite 
the  expansion  (4-32)  in  terms  of  the  local  Knudsen  number 


C  =  C 
D  o 


1+ 


H 


subject  to  the  condition 
~ -1 

K  A«l, 


(4-37) 


(4-38) 


then  lim  C±/CQ=  - (0.117±0.002)  becomes  independent  of  the  speed  ratio  S. 
Thus  the  apparent  increase  of  C±/C0  at  large  velocities  as  a  function  of 
S  is  a  consequence  of  the  fact  that  the  local  inverse  Knudsen  number  itself 
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becomes  proportional  to  the  speed  ratio  S.  Unfortunately,  the  condition 
(4-38)  implies  for  large  values  of  the  speed  ratio  that  K  ^«S  1  and,  hence, 
the  range  of  validity  of  the  expansion  decreases  with  increasing  values  of 
the  speed  ratio  S.  Of  course,  this  limitation  applies  to  all  theoretical 
results  obtained  by  a  Knudsen  number  iteration  procedure. 

Experimental  data  for  the  drag  coefficient  of  a  sphere  at  high  velocities 
have  been  reported  by  Kinslow  and  Potter  [40] .  In  terms  of  our  parameters 
these  data  correspond  to  SK  ^*=3  and  therefore  do  not  overlap  with  the  range 
SK  ^«1  of  our  calculated  values. 
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CHAPTER  V 

REMARKS 


In  this  report  we  have  considered  a  density  expansion  for  the  drag 
coefficient  in  terms  of  the  density  of  the  gas  stream.  When  the  mean 
free  path  of  the  molecules  is  sufficiently  large,  this  expansion  reduces 
to  an  expansion  in  terms  of  the  inverse  Knudsen  number.  Retaining  only  the 
first  two  terms  we  thus  obtained 


C  =  C  +  C  K 
D  o  1 


-1  , 


where  the  parameter  K  1  was  defined  as 


C5-1) 


K_1=\f2"  Trra2R  .  (5-2) 

The  coefficient  Cq  represents  the  drag  coefficient  in  the  limit  of  free 
molecular  flow.  It  was  shown  that  the  coefficient  is  determined  by  a 
set  of  well  defined  collision  integrals  associated  with  the  dynamical  motion 
of  two  molecules  in  the  presence  of  the  object.  These  collision  integrals 
can  be  formulated  for  objects  of  any  shape.  In  order  to  demonstrate  the 
feasibility  of  the  method  we  calculated  these  collision  integrals  for  the 
coefficient  of  a  disc  and  a  sphere  as  a  function  of  the  speed  ratio  S 
assuming  that  the  molecules  of  the  gas  are  reflected  diffusively  by  the 
object.  For  convenience  we  took  the  temperature  of  the  object  to  be  equal 
to  the  temperature  of  the  gas  stream  and  assumed  that  the  molecules  of  the 
gas  stream  could  be  treated  as  hard  spheres.  These  approximations  are  not 
essential.  Other  temperatures  of  the  object  and  molecules  with  more  com¬ 
plicated  interaction  potentials  can  be  handled  by  using  the  appropriate 


90 


AEDC-TR -74-79 


binary  collision  operators  T(1X)_. 

The  range  of  validity  of  (5-1)  is  determined  by  the  conditions 

K-1<<1  for  S^l  ,  (5-3a) 

and 

K-1«S-1  for  S^l  -  (5- 3b) 

In  order  to  extend  this  range  one  needs  to  consider  higher  order  terms  in 
the  .expansion  (5-1)  for  the  drag  coefficient.  In  analogy  to  the  density 
expansion  of  the  transport  properties  of  a  moderately  dense  gas,  one  may 
anticipate  an  expansion  of  the  form  [5,16,17,19)41] 

C=  C  +  C,K-1+  ClK-2lnif  1+  Cj2+  .  .  .  (5-4) 

o  1  2  2 

where  K  is  the  local  Knudsen  number  in  the  neighborhood  of  the  object. 

The  relationship  between  this  local  Knudsen  number  K  and  our  parameter  K 
was  discussed  in  Section  4.7.  In  view  of  (4-36b)  one  therefore  may 
anticipate  higher  order  terms  that  are  nonanalytic  in  terms  of  the 
parameter  K  1  as  well  as  in  terms  of  the  speed  ratio  or  the  Mach  number. 

For  objects  with  two  dimensional  geometry,  such  as  a  cylinder  or  a  strip 
whose  length  is  large  compared  to  the  mean  free  path,  these  nonanalyticities 
appear  already  in  the  first  inverse  Knudsen  number  correction.  However, 
the  nature  of  these  nonanalytic  terms,  as  well  as  their  practical  signi¬ 
ficance,  is  not  yet  well  understood  and  further  research  is  required. 
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*  APPENDIX 

EQUALITY  OF'  THE  CONTRIBUTIONS  PROM  RECOLLISIONS  AND  CYCLIC  COLLISIONS 

In  this  Appendix  we  prove  that  for  an  object  in  a  gas  stream  of  hard 
spherical  molecules  the  recollisions  and  cyclic  collisions  yield  identical 
contributions  to  the  drag  force,  i.e. 


EC1=ER1  ' 
EC2“  ER2  ' 


(A-la) 
(A- lb) 


where  ED,  and  E  are  given  by  (2-65)  and  E  and  E  by  (2-66) .  This 
identity  was  earlier  noted  by  Kuperman  in  the  zero  and  infinite  Mach  number 
limits  [19] ,  but  the  theorem  has  a  general  validity  independent  of  the  speed 
ratio  and  independent  of  the  detailed  nature  of  the  integrand.  It  thus  also 
applies  to  other  collision  integrals  pertaining  to  the  nearly  free  molecular 
■flow  regime,  such  as  those  derived  by  Kelly  and  Sengers  [18]  for  the  mass 
flux  to  a  droplet  in  a  supersaturated  vapor. 

In  order  to  prove  the  theorem  we  define  the  quantity 


g(v|’)=  jd A2  63(R1~R2+  VJTJ+  v^’t2)  Jdv|"n(v|'  ’|v”)  (v^'  '-v|') 


v|' -n2<o 


and  consider  the  integrals 


v1  •&  <o 
12  12  u 


V  /d»  12 


Vi2-d12<0 


(A- 2) 


IR^/d®12  ' 


(A- 3a) 


(A- 3b) 
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On  comparing  C2-65a)  with  C2-66a)  we  note  that  E =  E  . ,  if  I  =  I  . 

Cl  R1  R  C  • 

For  a  gas  of  hard  spheres  the  velocities  before  and  after  the  second 
collision  are  interrelated  by  (2-18) 


9 


(v*  'd 
1  12  12 


Let  us  introduce  a  new  unit  vector  6 '  defined  as 


&•  =  _  V12~ (V12,g12)ei2 


so  that 


v2,==  vi  “ 


(A-4a) 

(A-4b) 


(A-5) 


(A-6) 


We  choose  a  coordinate  system  with  the  Z-axis  in  the  direction  of  v\j^»  The 
polar  and  azimuthal  angles  of  in  this  coordinate  system  are  9  and  <|>  and 
vj2*ai2  =  cos6.  It  follows  from  (A-5)  that  the  polar  and  azimuthal  angles 
01  and  (J)1  of  o'  in  the  same  coordinate  system  are  0 '=—  ~9r  and  that 

A 

V12  "(T 1  =  cos0'=  -sin0.  Thus 


2TT  TT 

Ic=  jfdS;L2  ^V12‘&12Lg(v2,)=  v12  sin0  |  cos8  |  g(v^')  = 

^12*&12<0  °  */2 


=  v^2  j&$ '  j  de'sinQ1  |cos6'|  g(v2')  =  Jdfr'  S^v^').  (A- 7) 

0  %  v^2*d'r<o 
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On  comparing  (A-7)  and  (A- 6)  with  (A-3a)  and  (A-4a) ,  we  conclude  that  indeed 


(A— 8) 


“►  “►  4 

which  implies  the  equality  E  =  E  .  The  equality  E  =  E  0  follows  if 

R±  C  2,  R2 

4- 

is  identified  with  v  . 
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